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Abstract 

Multipartitional  extensions  of  Bell  (unipartitional)  polynomials 
are  shown  to  be  a  natural  and  strong  tool  in  the  study  of  multivariate 
compound  discrete  distributions  through  their  generating  functions. 
Modifications  of  exponential  polynomials  sin^ilify  proofs  in 
fluctuation  theory,  whereas  asymptotic  properties  of  such  polynomials 
are  used  to  establish  the  asymptotic  normality  of  a  wide  class  of 
combinatorial  distributions,  including  Stirling  and  C-numbers. 
Extensions  of  these  numbers,  the  non-central  Stirling  numbers  and  the 
mutli-parameter  Stirling  and  C-numbers  are  studied  in  conjunction  with 
distributional,  estimation  and  characterization  problems  related  to 
compound  distributions.  Combinatorial  and  occupancy-model  aspects 
are  also  discussed.  Diagnostic  tests  in  data  analysis  are  pointed  out. 

Keywords:  Multipartitional  Bell  polynomials,  non-central  Stirling 


numbers,  multiparameter  Stirling  and  C-numbers, 
combinatorial  distributions,  fluctuation  theory,  compound 
distributions,  characterizations,  regression  function. 


t 

OCCUPANCY  M2DELS,  BELL-TYPE  POLYNOMIALS  AND  NUMBERS 
AND  THEIR  APPLICATION  IN  PROBABILITY  THEORY 

1 .  GENERAL  BACKGROUND  OF  RESEARCH  AREA 

The  development  of  classical  probability  theory  has  its  origins  in  games  of 
chance  and  is  primarily  of  the  discrete  type.  Even  today,  it  is  cus ternary  to 
introduce  and  describe  many  discrete  probability  models  in  terms  of  urn  or 
occupancy  models.  The  models  themselves  and  their  ramifications  involve  concepts 
and  methods  of  combinatorial  analysis.  It  is  not,  then,  surprising  that 
combinatorics,  in  one  form  or  another,  has  always  been  in  the  foreground  of 
probabilistic  and  statistical-inferential  arguments  since  the  times  of  Pascal  and 
Fermat. 

However,  in  spite  of  an  abundance  Of  probability-related  combinatorial 
results  in  the  vast  mathematical  literature  over  a  long  period  of  two  or  three 
centuries,  many  of  the  results  are  "randomly"  and  sparsely  scattered  in  so  many 
books  and  professional  journals  so  that  they  either  remain  unknown  to  the  majority 
of  professionals  or  keep  being  rediscovered  all  the  time.  The  gap  was  widening 
also  due  to  the  fact  that  certain  closely  related  areas  of  mathematics,  such  as 
the  calculus  of  finite  differences,  were  neglected  in  modem  curricula.  It  is 
only  recently  that  attenpts  have  been  made  to  tie  up  such  results  in  a  systematic 
manner.  For  these  reasons,  monographs  such  as  Jordan  (1950) ,  David  and 
Barton  (1962)  and  Johnson  and  Kotz  (1977)  signify  a  good  move  towards  a  unified 
approach  of  the  methodology  in  the  wide-scope  area  of  ccmbinatorial  probability. 
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^  This  research  concerns:  occqpancy  or  urn  models.  Bell  or  parti  tier  polynomi¬ 

als  and  certain  kinds  of  numbers  (integers)  relating  to  Bell  polynomials  such  as 
Bell  lumbers.  One  can  arrive  at  these  numbers  fran  different  directions ,  as  will 
beccme  clear  in  the  sequel;  nevertheless  we  cover  all  such  numbers,  which  are  also 
connected  to  occupancy  models,  under  the  name  "Bell-type"  numbers. 

Occupancy  models  and  distributions  sure  well  understood.  Urn  models  con¬ 
stitute  the  usual,  though  not  the  only  •  .  way  to  introduce  most  discrete  distribu¬ 

tions.  The  recent  monograph  by  Johnson  and  Kotz  (1977)  provides  a  corrprehensi- 
ve  exposition  of  this  approach. 

Partition  polynomials ,  also  called  Bell  polynomials  by  Riordan  (1958  , 

1968)  are  connected  with  the  derivatives  of  a  composite  function  and  provide  a  power¬ 
ful  tool  in  the  treatment  of  combinatorial  and  probabilistic  problems.  Let 

A(t)  =  f(y(t)) 

and  define  =  d/dt,  Du  =  d/du,  (1) 

\  -  d"  A(t) ,  fn=  f  (u)  with  u=g(t) ,  eP  g(t)=  gR . 


where 


An  =  Yn  (f;  glf...,gn) 
y  ' s  a  partition  polynomial  defined  by 


YJf  .qi '  •  •  •  '9n^  *  ® 


n(n)  fe,!  •  -  .  fci! 


and  ike  summation  is  over  all  partitions  n(n)  of  n,  i.e.,  all  nen-negative  integers 
Jt| , . . . ,  kn  such  that 

+  2k^ . .  ,+nk^"  n  , 


k  -  +  kj  + 


represents  the  nunber  of  parts  In  a  given  parti ti^j.  The  So  called  exponential  po¬ 
lynomials/  En(g1 / . . .gn) ,  introduced  by  Bell  (1934  a,  b,),are  a  special  case  of  Y  , 


namely/  when  f.  ■  1,  k  ■  1,  2, ...  Thus 


“Yn(1;  91/*-*gn)  =  e"9  D"t  *  9=g(t)  (3) 

The  exponential  generating  function  (egf)  of  the  sequence  {v  n*o}(Y0=1) 
of  Bell  polynomials  can  be  written  in  the  form: 

or  n 

exp  (uY)  =  l  Y  (f;  g. , . . .g  ) 

n=0  n  '  n  ni 

k  ,  (4) 

00  u  g,  . 

=  exp  f(  — £T“)  =  exn^f  G  (u)J, 


where, in  the  exponential  expansions,  we  set 

G  (u)  =  exp[(ug)  -  gj  £%  ^  ,  gkS  gk,  Y n*  Yn.  (5) 

It  should  be  mentioned  that  (4) ,  in  conjunction  with  (5)  ,  attests  to  the 
general  fact  that  the  algebra  associated  with  egf ’s  is  what  is  known  as  the  Blis- 
sard  (or  symbolic  or  unbral)  calculus;  the  algebra  of  ordinary  generating  functi¬ 
ons  is  knoun  as  the  Cauchy  algebra. 

In  the  umhral  calculus,  a  sequence  a^  may  be  replaced  by  the  sequence 
a1  of  powers  and  when  all  operations  are  nerformed  the  exponents  are  changed 
back  to  indices.  For  f  ample,  if  A  (u) ,  B  (u)  and  C  (u)  are  the  egf ’s  of  the  se¬ 
quences  ,  b^  and  c^  ,  respectively,  and 

C  (u)  =  A  (u)  B  (u) , 


cn  m}  -OWk 
k=0 

so  that,  in  the  unhral  calculus  notation, 
cn  =  (a  +  b)n,  ar=an  ,  bn  =  bn  ; 


the  egf' s  behave  like  espcnential  functions; 


.  i 

U\ 

-  A  - 


C  (u)=  exp  (cu)=  exj^(a  +  b)=  exp  (au)  exp  (bu)=  A  (u)  B  (u) 

The  exponential  polynomial  E  is  closely  related  to  the  generat¬ 

ing  function  (gf) ,  Ri cardan  (1958), 


-.f.?  <k*'-V 


(6) 


of  the  number  C  kn)  of  nermutaticns  of  n  elements  with  k^  unit  cycles,  kj 

2-cycles,  etc.  is  the  cycle  indicator  of  the  syrrretric  group  and  is  expres¬ 
sed  in  terms  of  E  (t, , . . .  t  ) : 

n  i  n 

Cn(tV'tn)::En(t1'  fc2'  2*  t3'  *  *  *' (n-1) '  V* 

The  palynanials  are  connected  with  several  enumeration  problems.  For 
details, we  re*er  e.cr.,  to  Riordan  (1958,  1968)  .  We  mention  here,  e.g.,  that 
Cn  (t,...t)  is  equal  to  the  gf. of  the  signless  Stirling  numbers  of  the  first  kind 
c(n,k)  =  (-D^^stnjk) ,  where  s(n,k)  are  the  Striling  numbers  of  the  first  kind; 
c(n,k)  is  equal  to  the  number  of  permutations  of  n  elements  with  k  cycles;  similar¬ 
ly,  Cn(o,  t, ...,t)  is  the  gf  of  the  number  of  permutations  with  k  cycles  no  one  of 
which  is  a  unitary'  cycle.  Moreover,  the  polynomials  themselves ,  in  e- 

xactly  the  same  manner  as  the  C  ,  are  related  to  ordered  cycles  of  permutations. 

En  is  the  ordered-cycle  indicator  and  is  associated  with  Stirling  numbers  of  the 
second  kind,  S(n,k) .  Far  exarple,  the  polynomials 


(x)  =  E  (x, . .  .x)  =  E  S  (n,k)  xK 


n 


k=o 


(7) 


for  x=  1  give  the  Bell  numbers;  S^d )  equals  the  nunber  of  partitions  of  n. 

Another  application  of  En(g,j , .  ..g^)  with  g^=  (s)kx,  s>  o  or  an  integer, 
leads  to  the  polynomials 

Cn,s(x)=En(  (s)1x 


,...(s)  x  )  =  £  C  (n,  k,  s)  yk, 

n  k  =  0 
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where  the  C-nunbers  C  (n,k,s)  were  introduced  by  Cacoullos  and  Charalambides (1975) 
and  further  studied  by  Charalanbides  (1974,  1977) .  These  authors  shaded  in  a  se¬ 
ries  of  papers  (see  also  Cacoullos  1977)  that  Stirling  nunbers  of  the  first  and 
second  kind  and  C-nirnbers ,  as  well  as  certain  generalizations  of  these,  emerge, qui¬ 
te  naturally, in  a  unified  treatment,  via  egf's,  of  the  minimum  variance  unbiased 
estimation  problem  for  left-truncated  logarithmic  series,  Poisscn,  binomial  and 
negative  binomial  distributions. 

The  role  played  by  Stirling  numbers  in  occupancy  and  distribution  problems 
is  veil  known;  also  their  interpretation  in  terms  of  cycles  of  permutations  (Rior- 
dan  ,1 958).  Analogous  combinatorial  interpretations, along  with  their  probabilistic 
counterparts,  can  be  given  for  the  C-nunbers  and  their  generalizations.  If  s  is  a 
positive  integer  the  numbers C  (n,k,s)  are  associated  with  the  classical  coupons 
collector’s  problem.  If  s  is  a  negative  integer  then  n'.*C(n,k,-s)  /k'.  equals  the 
number  of  ways  n>  k  indistinguishable  balls  can  be  distributed  into  k  groups,  each 
of  s  cells,  so  that  each  group  contains  at  least  one  ball.  The  generalized  C-num- 
bers,  C(n,k,s,r) , correspond  to  the  situation  where  each  group  is  to  contain  at 
least  r  balls  (Charalambides,  1974) , 


The  above  type  of  nunbers  appear  in  other  contexts  as  well.  For  example, 
the  generalized  Stirlino  numbers  of  the  second  kind,  in  our  notation  S  (n,b,r) , coin¬ 
cide  with  the  numbers  B^’^for  j=o,  in  the  notation  of  Sobel  et  al  (1977).  They 

(n,r)  ft) 

were  introduced  in  the  study  of  the  Type  1-Diriclet  integral  1^  (r,n) .  This  is 

a  multivariate  extension  of  the  incarplete  beta  function,  1^,  far  multinomial-re¬ 
lated  probabilities .  In  fact) 

n  (b) 

S(n,b,r)  j£_  =  I^^n)  (8) 


I^^trjn)  is  wit  pr  ability  that  each  cell  in  a  binomial  distribution 


receives  at  least  r  balls.  Tables  of  S(n,b,r)  are  given  in  Sobel  (op.cit.)  for  b=1 
(1)  23,  n=b(1)  25,  r=1  (1)  [n/b]  .  Far  b=1  ,  the  incomplete  beta  function 

l5^(r.n)=I  (r,n-r+1)  =  ~- — — — -  f  xr1(1-x)n_rdx 
P  P  B (r,n-ti-r)  Jq 


can  be  written  as 


n-r 


I(1)  (r,n)=  Dr  l  ni±il  , 
p  j=0  r(r)jl 


whose  generalization  leads  to  a  Type  2 -Dirich let  integral  (see  Caccullos  and  Sobel, 
1966)  . 

Apparently,  the  first  systanatic  attempt  to  lose  Bell  polynomials  in 
probability  has  been  the  recent  study,  Charalanbides  (1977  )  of  compound (generali¬ 
zed)  discrete  distributions.  Let  N  and  be  independent  integer  -valued  randan  va¬ 
riables  with  probability  generating  functions  (pgf)  f  and  g,  respectively.  Then 
the  pgf,  P(t)  say,  of  the  compound  distribution  of 


Sn  X1  +  x2  +  •  •  •  +  Xjv 

is  given  by  the  compound  function  P (t)r.  f  (g(t))  and  ,in  view  of  (1)  and  (2)  and  the 
elementary  result 

pn  <f.9)  =  piVnJ  D"P(t)it=0  ' 


we  readily  obtain  the  compound  probability  function  in  the  form 


P0  (f,g) =  f(g(0)) 


pn  “Jr  yn(f!9l'-'9n) 


(9) 


where  fn=Du  1^^  •  Dn=Dtg^i  t=o. 


Fran  a  recurrence  relation  for  Y  one  obtains  the  basic  recurrence 

n 


pn+1 (f,q)=  n+T  if  gk,1f  Pn-k(f'g) 


fk=f. 


Similar  results  can  be  given  for  the  manents  of  S^* 


INTRODUCTION  AND  SUMMARY  OF  RESEARCH 


Within  the  general  scope  of  the  research  project,  the  results,  are 
exhibited  in  the  following  Sections  2-8.  Whenever  the  results  are  already 
published,  only  main  ideas  and  summaries  are  given  in  these  sections.  When¬ 
ever  the  work  is  still  in  the  process  of  publication,  it  is  presented  in  an 
appendix. 

Section  2  gives  some  general  theory  concerning  compound  (generalized) 
discrete  distributions.  It  also  discusses  ad  hoc  estimation  procedures  in 
conjunction  with  the  problem  of  modeling  certain  real  accident  data  by  using 
the  usual  discrete  distributions. 

Section  3  presents  extensions  of  Bell  polynomials  appropriate  for  the 
treatment  of  multivariate  compound  distributions. 

Section  1+  deals  with  the  asymptotic  normality  of  general  combinatorial 
distributions,  including  Stirling  and  C-numbers  as  special  cases. 

In  Section  5,  non-central  Stirling  numbers  of  the  first  and  second  kind 
are  defined  and  their  applications  in  convolutions  of  classical  discrete 
distributions  as  well  as  their  combinatorial  interpretations  in  occupancy 
models  are  discussed.  Another  extension,  the  multiparameter  Stirling  and 
C-numbers,  motivated  by  the  estimation  problem  for  multiply  truncated  power 
series  distributions  is  further  discussed  in  Appendix  A. 

Some  modifications  of  Bell  polynomials  useful  in  simplifying  proofs  in 
fluctuation  theory  are  given  in  Section  6. 

Finally,  Section  7  looks  at  compounding  from  the  point  of  view  of  mixtures 
of  distributions.  It  turns  out  that  the  regression  function  of  the  mixing 
variable  on  the  mixture  (compound)  variable  in  conjunction  with  identif iebility 
results  yields  characterizations  both  for  discrete  and  continuous  mixtures. 
Since  the  regressions  in  the  discrete  case  are  closely  related  to  Bell-type 
polynomials  and  numbers,  the  results  are  given  in  Appendix  B. 


-  8  - 


2.  MULTIVARIATE  DISCRETE  MODELS  GENERATED  BY  COMPOUNDING  ^ 

Multivariate  discrete  models  (N,Z)  generated  by  compounding  (generalizing) 

r\j 

an  integer-valued  positive  random  variable  (r.v.)  H  by  a  d-variate  discrete 
random  vector  £  =  (X^,...,Xg)  are  investigated  via  probability  generating 
functions  (p.g.f.)  and  it  is  shown  that  Bell  polynomials  and  related  numbers 
(e.g.  Stirling  and  C-numbers)  play  an  important  role,  not  only  in  expressing 
probabilities  and  moments, but  also  in  explicit  representations  of  the  conditional 
distribution  of  N  given  Z.  More  specifically,  the  (compound)  distribution  of 
Z  is  determined  by  the  representation 


Z 


(1) 


where  are  independent  observations  on  £,  which  is  assumed 

a1  a/ 

independent  of  the  generalized  r.v.  N. 

Several  applications  lead  to  such  models.  For  example,  N  may  represent 

the  number  of  car  accidents  in  a  given  locality  during  certain  time  period, 

X^  the  corresponding  injury  accidents,  X£  the  fatal  accidents,  X^  the 

injuries  and  the  fatalities.  The  joint  behavior  of  N  and  any  or  all 

of  the  X.  requires  an  (N,Z) -model. 

1  a, 

The  emergence  of  Bell  partition  (d=l)  and  multipartitional  (Section  3) 
polynomials  in  the  study  of  the  probabilistic  structure  of  an  (lv,Z)-model  is 
due  to  the  follwoing  basic  facts: 

1.  The  p.g.f.  G  of  (N,Z)  is  given  by 


G(u,v)  =  g. (u.g0(v) ) 

*\,  J-  C-  % 


(2) 


where  g^  is  the  p.g.f.  of  N  and  g2  the  p.g.f.  of  C;  the  p.g.f.  of 

the  compound  distribution  of  Z  is  the  compound  function  G,(go(v) )=G(l,v). 

%  1  d  n,  % 


2.  The  p.g.f.  of  N  given  Z  =  (Z^,...,^),  h2(u),  say,  ^s  given  by 


h  (u)  =  (u,0)/G^z)(l,0) 


where 


32og(u,v)  "1 
G  £  (a,f$)  =  - *■ 


^  32i...3Zk 

vi 


vk  1 


u=a,  za=z1  +  ‘--  +  2,  . 
’01  k 

Vi 


Some  general  results,  Cacoullos  and  Papageorgiou  [9],  concern  the 

conditional  distribution  of  M  given  Z=z,  when  d=l,  i.e.,  when  Z  is 

a  scalar  r.v.  It  is  shown  that  (for  details  we  refer  to  [9]): 

(i)  N|Z=z  is  a  convolution  of  n|z=0  and  another  nonnegative  r.v.  Y. 

(ii)  If  N  and  X^  have  power  series  distributions  (PSD),  then  the 

conditional  distribution  N|Z=0  is  the  same  distribution  of  N  with  a  new 

parameter  b^0  where  ^ is  the  parameter  for  N,  and 
r  R 

foOo)  -  L  is  the  series  function  for  X,. 

*  k^O  *  2  1 

(iii)  If  N  is  Poisson  and  X^  as  in  (ii),  then  the  r.v.  Y  of  has  a 

PSD  with  series  function  the  exponential  polynomial  ,2!bg8 , . . . ,z!bz0) , 

that  is,  a  general  combinatorial  type  distribution.  Harper  [21] . 

Fitting  some  actual  data  concerning  injury  accidents  (N)  and  fatal  accidents 

or  fatalities  (Z)  in  eastern  Virginia,  Leiter  and  Hamdan  [27]  used  a  Poisson- 

Bernoulli  and  a  Poisson-Poisson  model.  Negative  binomial -Bernoulli  or  negative 

2 

binomial -Poisson  models  gave  more  satisfactory  results,  as  Judged  by  the  x  - 
criterion,  Cacoullos  and  Papageorgiou  [lO] .  A  Poisson-binomial  model  [8],  was 
also  fitted  to  the  same  accident  data. 

In  addition  to  the  probabilistic  aspects  of  these  models,  inferential 
problems  are  also  examined.  Thus,  in  view  of  the  difficulty  in  obtaining 


explicit  solutions  of  relevant  maximum  likelihood  equations,  special  ad  hoc 


procedures  are  employed.  Such  are  the  methods  of  "even  points",  "zero  frequencies 
and  "ratio  of  frequencies". 

The  method  of  "even  points"  was  used,  [lo] ,  to  estimate  the  three 
parameters  (N,P  and  A)  in  the  NB-P  model;  in  addition  to  the  estimators 
X  and  Z  of  E(X)=NP  and  E(Z),  respectively,  use  is  made  of  the  equation 

G(l,l)+G(-l,-l)  =  2(Pee+P00) 
where  G( • » • )  is  the  p.g.f.  of  (X,Z)  and 

Pee  =  P[x=even,  Z=even]  ,  Pgg=p[x=odd,  Z=odd]  . 

This  yields  a  third  estimating  equation 

l+[Q+Pe"2X]"N  =  2(See+S00)/n  (Q  -  1+P) 

where  See  and  SQq  are  the  observed  frequencies  of  [x=even,  Z=evenJ  and 
[x=odd,  Z=odd]  in  a  sample  of  size  n,  respectively. 

The  method  of  "zero  frequencies"  uses  the  relative  frequency  fQO  in 
the  (0,0)  cell  and  the  proportions  f^  and  f  q  of  zero  observations  in 
the  two  marginals.  The  method  of"ratios  of  frequencies"  makes  use  of  ratios 
such  as  f^  q/ fflo  / fg  etc.  For  further  details,  we  refer  to  [9]  and  fio]. 

3.  MULTIPARTITIONAL  POLYNOMIALS 

Several  situations  (cf.  Section  2}  call  for  the  study  of  a  compound 
(generalized)  vector  random  variable  (r.v.) 

(1) 

where  X1,X„,...  are  independent  observations  on  the  r.v.  X  with  probability 
'v1  aA  'v 
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generating  function  (p.g.f.)  g(v),  say,  and  N  is  a  noifo-negative  integer- 

*  \ 

valued  r.v.  independent  of  the  X^,  with  p.g.f.  f(u).  Then  the  compound 
(generalized)  distribution  of  Y  has  p.g.f.  f(g(v)).  If  the  X.  are 

w  '*’1 

continuous  r.v.'s  with  common  moment  generating  function  (m.g.f.)  (or 
characteristic  function)  <p(v),  then  the  m.g.f.  (or  c.f.)  of  £  is  again  a 
composite  function,  f  (<p(v)).  Thus  we  could  consider  the  usual  properties 
of  %  (probability  function,  moments,  etc.)  in  terms  of  the  functions  f  and  g. 

The  case  of  scalar  r.v. 's  X^  led  very  naturally  to  the  use  of  Bell 
(partition)  polynomials  Yr  ,  Charalambides  [13J ,  since  Yr  may  be  regarded 
as  the  n-th  derivative  of  the  composite  function  A(v)=f(g(v))  in  terms  of 
the  derivatives  f^  of  f ( . )  and  g^  of  f ( . ) ,  k=l,2,...  .  Similar 
considerations  motivated  the  introduction  of  bipartitional  polynomials 
for  the  study  of  Y  in  the  bivariate  discrete  case,  Charalambides  [ ll* J . 
Essentially,  the  analogous  multipartitional  polynomials  can  be  used  for  the 
treatment  of  multivariate  distributions  of  (N,  Y)  (cf.  Cacoullos  and  Papa- 
georgiou  [9l  )  since  the  p.g.f.  of  (N,  Y)  is  f(u  g  (v))  as  stated  in 

'V  % 

Section  2. 

For  distributional  purposes,  it  is  convenient  to  define  a  bipartitional 
polynomial  in  terms  of  derivatives  as  follows  , 

Let  A(u,v)=f (g(u,v) )  and  set 


31  31 


V  ^  t-e(u,v)'  «u- 


m !  n !  f , 


m  n  g. 


Amn=Ymn=Ymn(f;g01’g10’,”,grnn)=^  kno!k,„!...k  !  (2) 


02  10  mn  i=0  J=0 


where  the  summation  extends  over  all  partitions  of  the  bipartite  indexes 

m  n 

(mn),  i.e.,  over  all  non-negative  integers  k. ,  satisfying  £  i  £  k.  =m, 
n  m  i=l  J*0 

.  I  J  I  k.  =n;  k  is  the  number  of  parts  in  the  partition.  The  expression  on 
’j=l  i-0  J 

the  F.H.S.  of  (2)  is  the  analogue  of  di  Bruno’s  formula  for  Bell  (partition) 


polynomials,  Riordan  [28]  ,  and  may  be  used^s  an  alternative  definition  of  Y^ 
(cf.  (2)  of  Section  l). 

An  immediate  consequence  of  the  fact  that  A  =Y  is  that  the 

mn  mn 

probability  function  of  Y=(Ylt  Y2),  with  p.g.f.  f(g(u,v)),  is 


Pm,n-PIil-nl,Y2“nl  *  Kin!  ’p01’p10 . . 


(3) 


s? f(t) 


‘-Poo’  Pij“PlXl=i »x2=j]>  poo'f(B(0'0))*r(poo> 


Moreover,  the  factorial  moments  \t,  X=E |(Y, ) - (Y~)  1  are 

Un,n;  *■  1  m  £  n* 


P(m,n)  Ymn(a;e(0,l)’e(l,0)*,**’e(m,n))  V“(k) 


(M 


with  and  6^  denoting  the  factorial  moments  of  N  and  (X^,X2), 

respectively. 

As  in  the  case  of  simple  Bell  (unipartitional)  polynomials,  using  the 

umbral  (Blissard)  calculus,  we  may  define  the  polynomials  Y  of  (2)  in 

mn 

terms  of  their  exponential  generating  function  : 

Y(u,v)  =  l  l  Y^  ^7- rr  =  exp[f{G(u,v)-G(0,0)}] 
m=0  n=0 

with  Afk,  G(u,v)  *  l  l  p-  TT  • 

*■  £_Q  J_Q  1  *  J  * 


This  can  be  used  to  derive  recurrence  relations  for  Y  .  Moreover  (5) 

mn 

(see  also  (k)  and  (5)  of  Section  1)  implies  the  following  important. 
Remark.  Bell  polynomials  can  be  used  equally  well  both  for  truncated  and 
non-truncated  versions  of  discrete  compound  distributions.  Truncation 
amounts  to  the  cancellation  of  certain  g.  ,  corresponding  to  truncated 

A  J 

values  of  the  r.v.’s. 

Bipartitional  polynomials  can  arise  in  another  situation  of  bivariate 


(5) 


compounding.  Let 


VXil+Xi2+'*'+XiN  i=1»2,  (6) 

1  1 

where.^  N^.Ng  are  non-negative  integer-valued  r.v.’s  with  joint  p.g.f. 

F(*»’),  the  X^j  are  independent  with  p.g.f. *s  g^(*)  (i*l,2)  and  the 

{X± j }  are  independent  of  (N^.Ng).  Then  the  p.g.f.  G  of  (Y^.Yg)  is 
easily  seen  to  be  the  composite  function 


G(u,v)  =  F(g1(u),g2(v)). 

The  special  case  fk=l  (k=0,l,..)  in  (5)  gives  the  bivariate  (bipartitional ) 
analogues  (gQpg.j_Q,g^,  *  *  ’  »gmn^  of  Bel*  exponential  polynomials 
(cf.  (2)  and  (3)  of  Section  1).  Thus  in  the  important  case  of  compound 
bivariate  Poisson  distributions,  when  either  N  or  (N-pNg)  is  Poisson, 
the  p.g.f.  of  (Y-pY,,)  takes  the  form 

G(u,v)  =  exp[x{g(u,v)-l}]  (7) 

in  the  former  case  (l),  whereas  in  the  latter  case  (6),  it  is  of  the  form 

G(u,v)=exp[xi(g1(u)-l}+X2{g2(v)-l}+X1^g1(u)g2(v)-l}]  (8) 

Two  examples  of  (7)  and  (8)  have  been  studied  by  Charalambides  and 

Papageorgiou  [l7] :  (a)  in  (l)  is  a  bivariate  binomial  i.e. 

g(u,v)=(p00+p10u+pQ1v+p11uv)n  and  (b)  the  Xp  in  (6)  are  independent 

ni 

binomials,  i.e.,  g^(u)  =  (p^u+q.^ )  .  They  provide  alternatives  to  Neyman 

type  A  models  used  by  Holgate  [22]  to  fit  certain  ecological  data. 

It  should  be  observed  that  the  exponential  polynomials  are 

associated  with  compound  bivariate  distributions  in  which  (Y^,Y2)  has  a 
p.g.f.  of  the  form 

G(u.v)  .  eh(u’T)  .  (9) 

In  this  situation  certain  useful  relations  may  be  stated. 

The  bipartitional  exponential  polynomials  ^mn^Soi’^10 


,...)  satisfy 


the  recurrence 


n  m 


Ea,n+1  *  I  I  •  E00=1  • 

*  6=0  r=0  * 

The  p. f.  P(m,n;h)=P[y^=m,Y2=n]  associated  with  (9)  is  given  by 


(cf  (3)) 


P(m,n;h)=e  '  Em#n^oi*h10*hll’  ’  * '  ,hmn^m:n: 


where  (see  (l)  of  Section  1) 


h  =  DBDr  h(u,v) j  _  n  . 
rs  v  u  ’  tu=0,  v=0 

As  regards  the  factorial  u(m,n;h)  of  (Y^.Yg)  we  have  (cf  ( U )) 


where 


p (m,n;h)  =  ^(c^.c^.c^ . *  *  * 


Crs  *  DvDu  h(u’v) lu=l,v=l  ' 


Using  (10)  we  get  the  recurrences 


P(m,n+l;h)  =  ~  l  1  "5^7^  P(m-r,n-s;h)  , 


n+1  Ln  ris! 

s=0  r=0 


p(m,n+l;h)  -11  (“)(")  cr  s+1P(®-r,n-s;h). 
s=0  r=0  ’ 


Marginal  and  conditional  p.f. 's  can  also  be  given  in  terms  of  E  .  Thus 


P^m]  =  eh(°’l)Em(h1 . hj/ml 


where  Em  is  the  corresponding  uniparti tional  polynomial  and 


hr  =  Djj  h(u,l)|u_0 

and  the  conditional  p.f.  of  Yg  given  Y^=m  is  given  by 


P  Yg=n  Y^=m  =  exp  h( 0,0)-h(0,l) 


Em,n(h01’h10>,,,,hmn 
n!Em^hl»"  •  ,hm^ 


For  more  details  we  refer  to  [ 14 ] 
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4.  ASYMPTOTICS  OF  COMBINATORIAL  DISTRIBUTIONS 

Relations  between  Stirling  and  C-numbers  and  Bell  polynomials,  as 
well  as  their  role  in  discrete-distribution  theory  and  occupancy. type 
problems  have  already  been  made  clear  in  the  preceding  sections. 

Another  interesting  probabilistic  aspect  of  Stirling  and  C-numbers 
is  tied  up  with  the  asymptotic  normality  of  the  so  called  combinatorial 
distributions.  A  first  result  in  this  direction  is  that  of  Harper  [2l] 
showing  the  asymptotic  normality  of  the  combinatorial  distribution  defined 
in  terms  of  Stirling  numbers  of  the  second  kind,  S(m,n),  namely,  the 
distribution  of 

PfXm=nJ  =  n=0,...,m,  m=0,l,...  (l) 

ID 

m 

where  B  =  T  S(m,n)  is  the  Bell  number.  Another  case  considered  by 

m 

n=0 

Charalambides . [l2] ,  is  when  S(m,n)  is  replaced  by  the  C-number 
C(m,n,s),  Cacoullos  &  Charalambides  [6].  In  a  general  combinatorial 
distribution,  the  S(m,n)  are  replaced  by  A(m,n)  which  are  assumed 
to  satisfy  the  "generalized  Pascal  triangle": 

A{m+l,n)=g(m,n)A(m,n)+h(m,n)A(m,n-l)  (2) 

where  g  and  h  are  positive;  also  the  r.h.s  of  (l)  is  multiplied 
by  An  ,  i.e.,the  general  combinatorial  distribution  is  defined  by 

P[X  =n]  =  \n  ,  A  (A)  =  l  A(m,n)An  .  (3) 

A  (A)  n=0 

XQ 

The  question  raised  here,  Kyriakousis  [26 j ,  is  under  what  conditions 
on  g(m,n)  and  h{m,n)  the  asymptotic  normality  (as  nr*®)  of  the 
general  combinatorial  distribution  obtains. 

In  particular,  the  problem  of  asymptotic  normality  is  studied 
through  the  asymptotic  behavior  of  the  polynomials 


ti 
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Am(X)  =  l  A(m,n)Xn  . 
n=0 

It  is  assumed  that  the  Am(X)  are  Bell  exponential  polynomials,  defined 
by  their  exponential  generating  function  : 

I  A  (X)  rr  =  exp{x[f(z)-f(0)j},  f(z)  =  l  c  z11. 

m=0  *  n£0 

For  example,  if  f(z)=ez,  then  A  (l)=B_,  the  so-called  Bell  number 

m  m 

and  A(m,n)=S(m,n) ;  if  f(z)=(l+z)S,  s>0,  then  A(m,n)=C(m,n,s) ,  the 

C-numbers.  Furthermore,  asymptotic  expressions  for  ratios  Am+jc( X)/Am(X) 

are  obtained  under  certain  conditions  on  g.  These  ratios  determine  the 

asymptotic  behavior  of  the  variance  (k=2)  and  the  mean  (k=l)  of  the 

X  .  Then  applying  the  normal  convergence  criterion  to  the  sequence 

{Xmn}  shows  the  asymptotic  normality  of  the  corresponding  combinatorial 

distribution.  Also,  a  useful  result  in  this  direction  is  that  of  Haigh  [ 20] , 

which  requires  to  show  that  the  polynomial  p.g.f.  of  X  has  real  roots. 

m 

Some  general  conditions  on  g(m,n)  and  h(m,n)  under  which  the 
corresponding  combinatorial  distributions  converge  to  the  normal  (as  nr*”) 
are  given  in  Kyriakousis  [26].  They  imply,  as  special  cases,  the  asymptotic 
normality  of  the  combinatorial  distributions  defined  by  : 
a)  the  signless  (absolute)  Stirling  numbers  of  the  first  kind  |s(m,n)|, 
which  satisfy  the  recurrence  (cf.  Appendix  A) 

|s(m+l,n) |=m|s(m,n) |+|s(m,n-l) | , 

(b)  the  binomial  coefficients  (®)  so  that  by  taking  X=p/q  in  (3)  one 

obtains  the  well-known  asymptotic  normality  of  the  binomial  distribution 

(c)  the  Eulerian  numbers 

E(m+l,n)=n  E(m,n)+(m-n+2)E(m,n-l) 

and  the  generalized  Eulerian  numbers,  Dwyer  [l9j ,E  (m,n),  a>0  integer,  with 
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*4 

recurrence 

S  E  (m+l,n)=(n+a)E  (m,n)+(m-n+2-a)E  (m,n)  , 

a  a  a 

(d)  the  Stirling  numbers  S(m,n)  of  the  second  kind  (a1(m)=bQ(m)=l, 
aQ(m)=0),  first  shown  by  Harper  |2l|  for  A=l;  also  the  C-numbers,  C(m,n,s), 
and  the  signless  numbers  jc(m,n,-s)|,  [12]  (see  also  Appendix  A).  Finally, 
the  non-central  Stirling  numbers  of  the  second  kind  S  (m,n)  ,  studied  by 
Koutras  [25],  with  the  recurrence  (see  Section  5  ) 

S  (m+l,n)  =  (n-a)S(m,n)+S  (m,n-l) , 

»  EL 

for  a<0  also  converge  to  normality 

5 .  NON-CENTRAL  STIRLING  NUMBERS  -  MULTIPARAMETER  STIRLING  AND 
C-NUMBERS. 

The  Stirling  numbers  of  th'  first  kind  s(n,k)  and  the  second  kind 

S(n,k)  are  usually  defined  t  ,  the  coefficients  in  the  expansion  of  the 

factorial  (x)n  in  powers  of  x  and  vice  versa  (see  also  Appendix  A). 

k  k 

The  non-central  ones  are  their  analogues  when  x  is  replaced  by  (x-c) 
for  some  real  c.  This  leads  to  an  equivalent  definition  in  terms  of 
exponential  generating  functions  (egf).  Thus  the  egf  for  the  non-central 
(around  c)  numbers  of  the  first  kind  sc(n,k)  is  found  to  be 

fkU)  =  l  8  (n,k)  ~  =  (l+t)C|ji-  log(l+t)|k;  (1) 

n=k 

similarly  the  egf  h(')  of  Sc(n,k)  is  given  by 

hk(t)  =  e"Ct  —■  (e^-l)^.  (2) 

These  expressions  can  be  used  to  obtain  relations  between  non¬ 
central  and  central  (c=0)  Stirling  numbers.  Moreover,  they  are  con- 


venient  for  the  definition  of  generalized  non-central  Stirling  numbers 
sc(n,k,r)  and  Sc(n,k,r) ,  by  subtracting  the  first  r  terms  in  the 
expansions  of  log  (l+t)  and  e\  respectively. 

Applications  in  probability  and  occupancy  problems.  Considering 
power  series  distributions,  (l-t)c,  c<0,  is  the  series  function  of  a  negative 

binomial  and  -log(l-t)  the  series  function  of  a  logarithmic  series 
distribution;  it  can  be  concluded  from  (l)  that  the  "signiless"  non¬ 
central  numbers  of  the  first  kind,  jsc(n,k)|  =  (-l)n-^sc(n,k)  are 
associated  with  the  convolution  of  a  negative  binomial  and  a  k-fold 
convolution  of  a  logarithmic  series  distribution;  similarly  the  generalized 
sc(n,k,r)  correspond  to  left-truncated  logarithmic  distributions. 

Another  use  in  probability  theory  of  sc(n,c)  is  in  expressing  the 

factorial  moments  x  of  a  random  variable  X  in  terms  of  its  moments 

n 

V  about  the  point  c,  that  is, 
x  ,c 

=  Jo  s=(n,klV‘=  • 

Thus  taking  c=y=E(X),  is  expressed  in  terms  of  the  central  moments 

Uy  .  Conversely,  the  Sc(n,k)  can  be  used  to  express  c  in  terms  of 
Tt n .  From  the  point  of  view  of  distribution  theory,  the  Sc(n,k)  are 
associated  with  convolutions  of  usual  and  zero-truncated  Poisson  distributions. 

An  occupancy-type  interpretation  of  Sc(n,k)  is  the  following  : 
n  distinguishable  balls  can  be  distributed  into  a  set  of  k  identical 
cells  and  ,}  distinguishable  boxes  so  that  every  cell  is  occupied  in 
S  j(n,k)  ways.  Similarly  for  Sc(n,k,r). 

For  additional  results  and  details  we  refer  to  Koutras  [25]. 


Multiparameter  Stirling  and  C-numbers. 


Another  extension  of  the  usual  Stirling  and  C-numbers  in  a  different 
direction  is  motivated  by  the  estimation  problems  when  several  independent 


samples  are  available  from  the  same  parent  distribution  but  the  truncation 
points  differ  from  sample  to  sample;  The  relevant  details  of  applications 
and  occupancy-type  interpretations  are  given  in  Appendix  A. 

6.  BELL  POLYNOMIALS  IN  FLUCTUATION  THEORY 


Bell  (exponential)  polynomials  have  been  used  in  the  study  of  generalized 
(compound)  discrete  distributions,  due  to  their  interpretation  as  derivatives 
of  composite  functions  (Sections  1  through  k) .  Another  important  area  of 
probability  theory  where  such  polynomials  provide  a  powerful  tool  is  the  so 
called  fluctuation  theory  as  developed  mainly  by  E.S.  Andersen  and  W.  Feller. 

Consider  a  generalized  random  walk  {S^},  where  Sjt*X1+*  •  *+Xk  k=l,2,... 
(Sq=0)  and  the  X^  are  i.i.d.  r.v.'s.  Using  a  classical  result  of 
Touchard  |32|on  the  number  of  permutations  of  n  elements  with  specified 
numbers  of  cycles  possessing  certain  properties,  in  conjunction  with  Spitzer's 
combinatorial  lemma,  yields  a  simple  combinatorial  proof  of  the  basic  result: 


Pn  =  PtSl  =  Q1,  Cn(&1,...,an) 

=  p[si<0,...,Sn<0]  =  Cn( 1-a^, . . . ,l-an) 


with  ak=P [sk>0] , 
indicator  function 
by  : 


k=l,2,...,n  (n=l,2,...)  and  Cn  denoting  the  cycle 

,  Riordan  [28],  related  to  Bell  polynomials  Yn(t^,. 


cn(tr. 


'  =  Yn*tl,t2,2!t 


.3,., 


(n~l)!tn). 


Further  simplifications  are  obtained  by  introducing  the  polynomials. 


Ck,n(x»y)=Ck,n(xl’* 


’xn’yl* ' 


,y„)=0  C.(x)c, 


n-k 


(y). 


Exploiting  certain  properties  of  n(x,y),  one  can  show  the  basic  result 
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and  hence,  e.g.,  the  recurrence  ^ 

n-k 

plN„4.i=kJ  =  — T7T  l  U-a  ,  i )  P[N  =k]  , 

A  n+1  ■>  n-l+k  Ln  r+1  1  n-r  J 

r=0 

N  denoting  the  number  of  positive  partial  sums  S,  ,  k=l,...,n. 

n  k 

Also  simplified  proofs  can  be  obtained  of  some  results  concerning 
symmetrically  dependent  (exchangeable)  r.v. 's  after  proving  a  result 
along  the  lines  of  Spitzer’s  combinatorial  lemma.  For  details  we 
refer  to  [l6] . 

7.  CHARACTERIZATIONS  OF  COMPOUND  DISTRIBUTIONS  BY  REGRESSION 
AND  BELL-TYPE  POLYNOMIALS  AND  NUMBERS. 

So  far  Bell-type  polynomials  were  discussed  in  relation  to  the 
distribution  of  a  generalized  (compound)  discrete  random  variable  (rv); 
also  in  generalized  random  walks  (fluctuation  theory).  In  the  latter 
case,  it  was  indicated  how  Bell  polynomials  can  be  extended  to  provide 
simple  proofs  for  some  known  basic  results.  In  the  former  case,  the 
probability  generating  function  of  the  generalized  rv  Y=XvZ, 

Y  =  Z1+Z2+"'+ZX’  (1) 

can  be  expressed  in  terms  of  Bell- type  polynomials.  For  example, 
Stirling  type  polynomials  (with  coefficients  Stirling  numbers  of  the 
second  kind)  appear  whenever  the  generalizing  variable  Z  is  a  Poisson 
and  C-type  polynomials  (with  coefficients  C-numbers)  whenever  the 
Z-variable  is  either  a  binomial  or  a  negative  binomial,  [9]. 

Bell-type  polynomials  can  also  be  used  to  express  the  regression 
(posterior  mean)  m(y)  =  E [ X ] Y  =y ]  of  a  compound  mixing  r.v.  X  on  the 
mixture  variable  Y,  with  (absolute)  probability  function  (p.f.) 

p(y)  =  I  p(yjx)  f(x) 
x 


(2) 
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(p( •  j x )  denoting  the  conditional  p.f.  of  Y(X)  ,^given  X=x  and  f ( • ) 
the  prior  p.f.  of  X).  This  came  up  in  specific  applications  of  the 
following  general  characterization  result. 

Let  p(y|x)  in  (2)  be  a  binomial,  negative  binomial  or  Poisson 
for  every  x=l,2,...  .  Then  m(y)  characterizes  both  p( • )  and  f ( ' ) . 
(The  details  are  given  in  Appendix  B, exhibiting  also  several  examples 
of  Bell-type  polynomials  and  numbers). 

In  fact,  m(0)  alone  is  sufficient  to  characterize  p(*)  and  f(*) 
provided  p(*  jx)  is  an  x-fold  convolution  of  the  same  r.v.  Z  for  every 
s  =  p[Z=o]>0,  Cacoullos  [5].  Clearly,  by  varying  arbitrarily  the  X- 
distribution,  an  infinite  variety  of  bivariate  discrete  distributions 
can  be  characterized  in  this  fashion.  The  result  is  based  on  the  fact 
that  the  p.g.f.  h  of  X  satisfies  the  differential  equation 

m(0)  =  m(0;s)  =  sh'(s)/h(s)  0<s<l  . 

Hence  f ( • )  is  determined  and  by  (2)  also  p('). 

Finally,  characterizations  in  terms  of  m(y)  are  obtained  for 
continuous  analogues  of  (2),  when  X  and  Y  are  continuous.  Here  the 
role  of  X  is  played  by  a  continuous  parameter  0,  which  is  itself  a 
r.v.  with  some  prior  density  f { • ) .  However,  since  these  results  fall 
rather  outside  the  main  scope  of  this  research,  no  further  details  are 
given  here. 
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ABSTRACT 

Multi parameter  Stirling  end  C-numbers  are  defined  via  exponential 
generating  functions  and  basic  recurrence  relations  are  given;  also,  some 
combinatorial  and  occupancy  type  interpretations  are  provided. 

Recurrence  relations  are  derived  for  certain  ratios  of  simple, 
generalized  and  mul ti parameter  Stirling  and  C-nuribers.  Tnese  recurrences 
are  useful  in  the  computation  of  minimum  variance  unbiased  estimates  (mvue) 
for  classical  discrete  distributions  truncated  on  the  left.  Asymptotic 
relations  betveen  these  numbers  are  also  included. 


Key  words:  multi parameter  Stirling  and  C-numbers,  exponential  generating 
functions,  recurrence  relations,  mvue,  left  truncation,  power  series  distributions. 
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1 .  INTRODUCTION 

i 

The  Stirling  numbers  of  the  first  and  second  kind  are  .jess  known  among 

statisticians  than  among  people  dealing  with  combinatorics  or  finite  differences. 

Only  recently  have  they  made  their  appearance  in  distribution  theory  and  statistics. 

They  emerge  in  the  distribution  of  a  sum  of  zero-truncated  classical  discrete 

distributions:  those  of  the  second  kind,  S(m,n),  in  the  case  of  a  Poisson  distribution 

truncated  away  from  zero,  Tate  and  Goen  (1958),  Caeoullos  (I96l);the  signless 

(absolute-value)  Stirling  numbers  of  the  first  kind,  |s(m,n)|,  in  the  logaritlimic 

series  distribution,  Patil  (1963).  In  general,  such  distributional  problems  are 

essential  in  the  construction  of  minimum  variance  unbiased  estimators  (mvue)  for 

parametric  functions  of  a  left-truncated  power  series  distribution  (PSD). 

Analogous  considerations  for  binomial  and  negative  binomial  distributions 

truncated  away  from  zero  motivated  the  introduction  of  a  new  kind  of  numbers,  called 

C-numbers  by  Caeoullos  and  Cbaralambi des  (1975)-  These  three-parameter  C-numbers, 

C(m,n,k),  were  further  studied  by  Cbaralambi des  (1977).  "who  gave  the  representation 
a 

C(m,n,k)  =  ][  k  s(m,r)  S(r,n)  in  terms  of  Stirling  numbers  of  the  first  kind, 

r=n 

s(m,r),  and  the  second  kind  S(r,n).  Interestingly  enough,  this  representation  in 
a  disguised  form  was,  in  effect,  used  by  Shumway  and  Gurland  (i960)  to  tabulate 
C-numbers,  involved  in  the  calculation  of  Poisson-binomial  probabilities. 

The  so-called  generalized  Stirling  and  C-numbers  emerged  as  a  natural  extension 
of  the  correspondi ng  simple  ones  in  the  study  of  the  mvue  problem  for  a  PSD  truncated 
on  the  left  at  an  arbitrary  (known  or  unknown)  point  Charal abni des  (I97^b).  It  should 
l>e  mentioned  that,  in  particular,  the  generalized  Stirling  numbers  of  the  second 
kind  were  independently  rediscovered  and  tabulated  by  Sobel  et.  al  .  (1977),  in 
connection  with  the  Incomplete  Type  I  -  Dirichlet  integral. 

Tlie  mult  iparamet »  r  Stirling  and  C-numbers  are  the  anBlogu'-s  of  g< n<  rnl  i  st  c 


Stirling  and  C-numbers  in  a  multi-sample  situation  where  the  underlying  PSP  is 
multiply  truncated  on  the  left,  Caeoullos  (1975)*  (1977). 
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Recurrence  relatione  for  ratios  of  Stirling  and  C-numbers  are  necessary^ 
because  the  mvue  of  certain  parametric  functions  of  3 eft-truncated  logarithmic 
series,  Poisson,  binomial  and  negative  binomial  distributions  are  expressed  in 
terms  of  such  ratios.  These  recurrences  bypass  the  computational  difficulties 
which  come  from  the  fact  that  the  numbers  themselves  (but  not  the  ratios  of  interest) 
grow  very  fust  with  increasing  arguments.  Recurrences  for  ratios  of  simple  Stirling 
numbers  of  the  second  kind  were  developed  by  Berg  (1975). 

The  main  purpose  of  this  paper  is  to  provide  recurrences  for  certain  ratios 
of  multi  parameter  Stirling  and  C-numbers,  thus  unifying  several  special  results, 
includir'7  those  of  Berg  (1975).  For  the  development  of  the  topic,  we  found  the 
use  of  exponential  generating  functions  (egf)  most  appropriate,  both  for  introducing 
the  numbers  themselves  and  deriving  recurrences  as  veil.  Without  claiming  completeness, 
we  included  certain  basic  recurrences.  el ^t-rvea  elsewhere,  Cacoullcr  (1975)  (1977), 
it  is  emphasised  here,  once  more, that  in  the  study  of  PSD's  the  egf  approach  is  the 
one  suggested  by  the  probability  function  itself  in  its  truncated  form.  Also,  we 
found  it  appropriate  to  include  certain  asymptotic  relations  between  Stirling  and 
C-numbers,  which  reflect  corresponding  relations  between  binomial  and  Poisson 
distributions  or  logarithmic  series  and  negative  binomial  distributions. 

A  typical  result,  which  involves  ratios  considered  here,  is  the  following. 


Let  x..,  t=l,...,n.  he  a  random  r.amnle  from  a  left-truncated  one-parameter  PSD 

l  J  l 


distribution  with  p.f. 

a. (x)6X 

p(x;6)  =  f.(e,r'J  •  x=ri’ri+1”*- 

CO 

where  f.(6,r.)  =  £  a.(x)6X,  i=l,...,k.  If  the  truncation  point 

1  x=r. 

i 

r=(r, , . . .  ,r.  )  is  known  and  a.(x)>0  for  every  x>r .  ,  i=3,...,k,  then, 

i  K  J  -1 


d.i) 


Cacoul 1  os  (1977),  for  every  J=l,2,...,  0j  is  estimable  and  its  (unique)  mvue, 
hared  on  all  k  independent  samples  {x^},  is  given  by 


* 


-Ah  - 


a{n.-J;n,r) 

e  (m)  =  (m)  -7 - v — 

J  J  a(m;n,r) 


where-  n=(nlt . . .  .n^) ,  r=(r1, .  . .  ,rR) ,  (m)j=m(ni-l)  •  •  •  (m-J+l) 


and 


^(1.2) 


l(ir‘U2 »£>  =  rr~T  I  D  f1  ^(x  ), 
nl*'"V  i=l  .1=1  1 


where  the  summation  extends  over  all  ordered  K-tuples  (li=n^+’  •  ,+nj,  ) 


of  integers 


(1.3) 


k  n.  - —  -  nati cfying 

x.  ->.r.  ,  £  ]>  x.,=m.  In  the  cases  of  interest  (Poisson, binomial ,  etc), 

J  i=l  J=1  J 

the  numbers  (integers)  a(m;u,r)  turn  out  to  be  Stirling  or  C-numbers,  depending 


on 


the  series  function  f^  in  (l.l),  which  at  the  same  time,  suggests  the 
corresponding  egf  of  these  numbers. 


2.  MULTI ^PARAMETER  STIRLING  NUMBERS  OF  THE  FIRST  KIND;  DEFINITION- 
GENERAL  PRO P E  RT 1 ES . 

Let  r1,...,rJf  and  n^,...,^  be  non-negative  integers  (k>l).  The 
mul 1 i parameter  Stirling  numbers  of  the  first  kind  with  parameters  r=(r,  ,rot . .  .  ,r.  ) 
and  ^=(n1,n2»...,nJr),  to  be  denoted  by  s(m;n,r),  can  be  defined  (cf.  Cacoullos,197^) 
by  the  egf 


=  I  s(m»s m :  =  I  T  ~r 

~  m=r 'n  i=l  i * 


log(l+t ) 


-  (-D-5-1  f 

J=1  JJ 


(2.1) 


where  ve  set  m=r'r)  =  r^n^-i-  *  *  •♦r^iij ,. 

The  special  case  k=l,  r^=r,  n^=n  yields  the  generalized  Stirling  numbers  of 
the  first  kind,  s(m ;n,r),  defined  by  Charal ambi des  ( 197^ a)  ,vhi 1 e  k  =  l,  r=l 
gives  the  simple  Stirling  numbers  of  the  first  kind,  s(m,n).  Froposi t ions 
2. 1-2. 3  summarize  basic  properties  and  recurrences  for  r(ir;n,r)  hnd 
facilitate  their  computation . 

Remark  2._1_.  In  Die  sequel,  in  order  to  avoid  unnecessary  compl  i  rati  on:  is  t:.i 
recur  rt  neer ,  vt  assume  that  all  n^>C';  some-  n,.  ,  say  v,  are  zero,  then  t!.(  j  ammeter 
k  becomes  k '-k-v  and  t  he  necessary  m  »di  fi cat  ions  art  obvious. 
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Proposition  2.1.  The  multi  parameter  Stirling  numbers  of 
the  first  kind  s(mjn,r)  have  the  following  representation 


n  • 

k  i 


s(m;n,r)  =  (-1) 


rnn  f 

1  km  i=l  3=1  13 


(2.2) 


where  N  =  n^+...+n^  and  the  summation  extends  over  all  ordered  K-tuples 

of  integers  x.  satisfying  the  relations  _ 

"  k  ni 


v..  >  r.  ,  i  =  1 . k  and 

13  —  1 


Proof.  We  have 


tt'V 

i=l  3=1  3 


Y(t  ,r^ )  =  log(l+t) 


r.-j 

•£ 


v  1  ***  vi  +  ^ 

(-1)  1  ~  =  JT  (-D  .  i=l . k.  (2.3) 


Forming  the  Cauchy  product  of  series,  we  find,  in  virtue  of  (2.1), 

«.<'*>  fh*  =n  [1*,^-  t  fi^r  (*.« 

-  i=l  i=l  L  J  m=r'n  m  i=l  3=1  *13 

where  ^  has  the  same  meaning  as  abover 


Comparing  (2.4)  with  (2.1)  we  get  (2.2). 

To  obtain  recurrence  relations,  we  make  use  of  the  easily  verified 

difference-differential  equation,  satisfied  by  the  egf  g  (t,r)  in  (2.1), 

n  ^ 


namely. 


a«>  jj  en(t;r)  =  t  tvl  »s-.1“*r>  ,2-5> 

-  i=l 

where  e^  =  (0,..., 0,1,0 . 0),  i.e.,  a  k-com.ponent  vector  with  zero 

components  except  the  i-th  component  which  is  equal  to  one. 

Proposition  2.2.:  (m,n)  -  wise  relations:  The  numbers  s(m,n,r)  satisfy  the 
recurrence  relation 

k  r.  -1 

s(m+ltn,r)+m  s(m;n,r)  =  Y~  (-1)  1  (m)  ,  s(m-r.+l*,n-e.  ,r)  (2.6) 

-  -  **  ~  r — -  r.-j  1  -  “1 

1=1  1 

with  initial  conditions 

s(0;0,r)  =  1,  s(0;n,r)  =  0  whenever  )~r.n.>  0,  s(m;n,r)  =  0  if  m  <  r'n. 

*  ^  ^  f  \  lx  •  *■' 


Troof .  Fquation  (2.5)  in  virtue  of  (2.1)  can  be  written  as 


mtr.  -1 


JL.  m-1  •  i.  •_  r .  -1  i 

)  s(m;n,r)  =  L  L,  s(m;n-e .  ,r)  f - 

m-.r’n  i  =  l  n,=r’n-r. 


(2.7) 


Equating  the  coef f i cients  of  t  ‘/m!  in  (2.7)  yields  (2.6).  Note  that 
equation  (2.6)  for  k  =  1,  r^=  1  ^ives  the  well-known  recurrence  for  the 
simple  Stirling  numbers  of  the  first  kind 


s(m+l,n)  =  s(m,n-l)  -  in  s(ir,,n). 


(2.6) 


Proposition  2.3.  (m;n,r)  -  wise  relations:  The  numbers  s(m.;n,r)  satisfy 

n±~  '  jr.  (m)Jr 

s(m;n,r+e.)  =  ) (-1)  - 5 — t—  s(m-  jr .  Jn- je.  ,r)  ,  i=l,...k  (2.9) 

1  j=0  j!  (r.)3  1  ~  -1  ~ 

Proof .  We  have,  using  also  (2.3) 


r .  n .  ,  i 

=  d^[T<,;-r),<-1)  -f— ]  n 

-  i  u  l  3=1  3  J  J 


(2.10) 


and  using  the  binomial  expansion 
r,-  _  n,-  n- 


r.  a-i  “i  111  /n.\  Dr.  J  i 

[V(t  ,ri)  +  (-l)  2  =  Z-  {,  J(n.-j):  gn  _.(t,r.)(-l)  11_^_(2.11) 

L  i  J  j  =  0  J  j  J  r. 


we  can  write  (2.10)  as 


.. .  j.  .  _ —  m.+  jr . 

t,  h  =  £ -.7”f — 2_,  .  E(n'^i!i’r) '-n::—  (2-22) 

n.=r'n+n.  3  =  0  3;  r.  r;=.r'n-3r. 

-~i  1  -  ~  1 

hence,  equating  the  coefficients  of  t^'/m!  we  obtain  (2.9). 

Signless  rnultiparameter  Stirling  numbers . Trom  the  recurrence  relation  (2.6), 
it  follows  that  the  numbers  s(m;n,r)  are  integers.  Moreover  from  the 
representation  in  (2.2),  we  conclude  that  s(m;n,r)  is  an  integer  with  sign 
(-1)*'  r‘,  where  N  =  r.^+...+n^.  Therefore  if  we  multiply  (2.6)  by  (-l)r'  +  1  , 


we  obtain 


I  s(r,+  l  ;n,r)  I  =  ir.j  s(r.;n  ,r)|  |  s(n.-r.+l*,n-e^  ,r)  |  (2.13) 

we  call  js(m;n,r)|  the  signless  (positive)  multiparamieter  (k-parameter) 

Stirling  number  of  the  first  kind.  We  will  show 

Proposition  2.4*  The  egf  of  |s(m,n,r)|  is  given  by 

k  ~  r-“3  •  _ n . 

f^v»r)  =  |  s(ir.;n,r)  |  ~  =  T|  “T  [-logd-t)-  J_  4-  J  (2.1<0 

-  r  =r’n  ’  i=l  i"  L  j  =  l  J 

Proof.  Prom,  the  difference  equation  (2.13)  it  is  easily  verified  that  the 
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egf  £*(tjr)  satisfies  the  difference-differential  equation 

(1-t)dT  hn(t;-r)  *  jlj  1  ' 
vhich  in  turn  yields  (2.ll»). 


(2.15) 


Alternatively,  (2.lk)  leads  to  the  representation  of  |s(m;n,r)|  as 
obtained  from  (2.2). 

3.  RATIOS  OF  MULT I PARAMETER  STIRLING  NUMBERS  OF  THE  FIRST  KIND. 


We  define,  as  ratio  of  multi  parameter  Stirling  numbers  of  the  first 
kind  vith  respect  to  argument  m,  the  function 


s(m+l;n,_r) 

R1(m;n,r)  =  — 

1  -  s(m;n,r) 


(3.1) 


Ratios  vith  respect  to  the  arguments  n^,^,  i  =  1,  —  ,k  can  also 
be  defined.  The  main  reason  for  considering  ratios  vith  respect  to  m 
is  seen  from  (l.l),  vhich,  actually,  involves  reciprocals  of  R^,  vhen 
ve  are  concerned  vith  the  parameter  of  a  logarithmi c  series  distribution. 

Proposition  3.1.  A  recurrence  relation  for  the  ratio  R^(m;n,r), 
independent  of  the  multi  parameter  Stirling  numbers  of  the  first  kind, 
is  given  by 


.  (m)  .  r.n,  ..  , 

k  r  -1  J  J  r-j-xji 


I  - /-"Vt - t  I  R.  (m-r  ,+l-i ',n-e  ,r) 

'T  n'ri  1  •) 

R^(ir.;n,r)+m  =  x - (3.2) 


1  T  R  (m-i;n,r) 
1=1 


-  A8  - 


for  n  ►  1  and  ir,  >  r'n  ,  with  Die  boundary  conditions 


Rj(m,l,r)  =  -ir. 

and  , 

k  n.  fc  n 

n-1  I  —  x  — 

1~2  -3'1  r  J  (r.+l)  n  r  1 

J  J  i  =  1 


(3.3) 


(3.h) 


Proof .  Using  equation  (3.1),  it  can  be  easily  seen  that 
tri -r'n  s(m;n,r) 


ri-X'fi  s(m;n,r) 

IT  R,(m-i,*n,r)  =  — r~~ — \ —  • 

'  '  l  s(r’n,n,r) 

1  1  •>  ^  ^ 


(3.6) 


But  equation  (2.2),  for  m  =  r'n,  m. ,  =  m.-  =  . . .  =  ti.  =  r.,  becomes 

^  in .  1 


s(r’n,n,r)  =  (-1)-  -  K 


(r'n).' 


k  Jk  n. 

n  v  n 


i=l  i=l 


Consequently,  equation  (3.5)  becomes 


s(it.;n,r)  = 


n.r  n 

l'H~K  fr'nV  TT  R  (m-i;n,r) 
r^i )  lr  n}-  j=i  1 _ 

&  •*’  IS  ‘ 


Tror.  equations  (2.6)  and  (3.1)  we  have 

k  r.-l 


R  (n.;n,r)+ir.= 


£~(-l)  3  (m)  s(m-r.+i;n-e.  ,r 
: — -  r.-l  1  ~  — 
3=1  3  J  J 


s(m;n,r) 


(3.6) 


(3.7) 


(3.6) 


and  substituting  for  sfr.-r^  +  l  ,n-e_.  ,r)  and  s(m,n,r)  from  (3.7)  yields  (3.2) 


By  definition 


s(r 'n+1 ,n,r) 

Y^*^)=  'ITP^r/r  ’ 

using  equation  (2.2)  for  m  =  r’n  +  1  ,  =  n^j, 

:‘iVV  for  1  1  ] . V 

required  formula  (3.^)  is  easily  obtained. 

The  special  case  k  =  1  yields 


(3-9) 


1  •  '  “  JT*#  —  UJ ,  «  -"*•••  “ 

1,1-1  1,1+1 

and  equation  (3.6),  the 
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Proposition  3.2.  A  recurrence  relation  for  the  ratio  R3(m,n,r),  in¬ 
dependent  of  the  generalized  Stirling  numbers  of  the  first  kind,  is  given  by 


rn  (m) 


r-1 


m+l-rn 


R3(m,n,r)+m  = 


(rn). 


7T  v 

»=i  _ 


m-r+l-i ,n-l ,r) 


m-rn 

]*]'  R1(m-i  ,n,r) 
i=l 


(3.10) 


for  n  £  1  and  m  >  rn,  with 

rn(rn+l) 


R^(m,l,r)  =  -m 


RJl(rn,n,r)  =  -  r+3 


(3.11) 

(3.12) 


Also  for  k  =  1,  r  =  1  we  obtain,  ^  -  t 

Proposition  3.3.  A  recurrence  relation  for  the  ratio  R^trn.n),  independent 
of  the  simple  Stirling  numbers  of  the  first  kind,  is  given  by 


R3(m,n)+m  = 


iri+l-n 

TT  R  (m-i,n-l) 


i=l 


(3.13) 


Rjlm.n.r)™  = 


R3(m-l,n,r) 


for  n  >  1  and  m  >  rn.  R3(m,l,r)  and  R3(rn,n,r)  are  given  by  (3.1i)  and 
(3.12),  respectively. 


Proof .  Using  equation  (2.6)  with  k  =1,  we  have 


R3(m,n,r)+m  = 


(-l)r  1  (m)  s(m-r+l,n-l,r) 
r- 1 


s(m,n,r) 

from  which  equation  (3.16)  can  be  easily  derived. 


(3.17) 


Applying  Proposition  3.4  with  r=l  gives 
Proposition  3.5.  An  alternative  recurrence  relation  for  the  ratio  R^tm.n) 
is  given  by 


R3(m,n)+m 


[R3(m-l,n)*m-l]  R^ (m- 1 ,n-l ) 


R3(m-1 ,n) 


n>l ,  m  ">  n . 


m-n 

][  R.(m-i.n) 

i=a 

for  n  >  1  and  m  >  n,  with 

R3(m,l)  =  -m  (3.14) 

R^(n,n)  =  -  n(n+l)/2  (3.16) 

Proposition  3.4.  An  alternative  recurrence  relation  for  the  ratio  R3(ir.,n,r) 
is  given  by 

jR  (m-l,n,r)+m-lj  R  (m-r,n-l,r) 
m  1  1  1  -  (3.16) 


(3.16) 
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'4.  MULTI  PARAMETER  STIRLING  NUMBERS  OF  THL  SECOND  KIND 

The  mult iparameter  Stirling  numbers  of  the  second  kind  S(ir.;n,r) 


defined  by  their  egf 


r.-l  .  .  n. 


fn(t;r)  =  r  s(n;n.r)  ^  J  ^  f 

-  m^r'n  i-l  a  j=0  J  J 


( L  .  1 ) 


Taking  k=l,rT=r  gives  the  generalized  Stirling  numbers  of  the  second  kind,S(n,n,r 
Charalambides  (l97^a),  taking  k=l,  r=l  defines  the  simple  Stirling  numbers 
S(m,n).  Tlie  following  properties  of  S(ir.,n,r)  can  easily  be  established  (of. 
Section  2). 

a)  They  have  the  representation 

k  n ' 

S(m;n,r)  s  H  n  n 

V***V  m  i=l  j=l  xi y 

where  the  summation  extends  over  all  ordered  N-tuples  (i.^n,4 •  •  -+nj,)  of 


integers  x- .  satisfying 

•*  J  ,  n . 

k  i 

x..  >  r.  ,  i  =  l,...,k  and  Y~  \  x.  .  =  m. 
lj  —  l  ^  4-^  i j 

b)  They  satisfy  the  following  recurrence  relations 


S(n+i;n,r)  =  N  S(m.;n,r)  ^  T  /”  )  S(m-r  .tl ;n-e .  ,r) 

—  -  -v  frj  \  r^-l  J  l  -l  - 

n .  (m) . 

C2-  n  3ri 

S(m;n ,r+e .  )  =  2 _  (-T)J  — - r  S(m- jr.  ;n-je.  ,r) 

j=0  j!(r.!)3  i- 


(U.3) 


with  initial  conditions 


S(0;0,r)  ~  1,  S(0;n,r)  =  0  whenever  J  r.n.>0  and  S(m;n,r)-0  if  m<r 'n . .  5) 


These  follow  from  the  difference-differential  equation 

3t  Vt!£>  *  *  fD(t>r)  *  t  ;r)  tri  /irj-D! 

-  ~  1~1  ~  -l 


(■4.6) 


It  can  be  easily  seen  that  the  representation  (*4.2)  prr  ’des  the  following 

comhinat orial  interpretation  in  terms  of  occupancy  numbers. 

Proposition  4.1.  The  number  of  ways  of  placing  m  distinguishable  balls  into 

N  =  n,t...tn,  cells  so  that  each  cell  of  the  i-th  group  of  n.  cells  contains 
J  k  c  ‘  i 

at  least  r.  balls  for  i  =  l,...,k  is  equal  to  n  J . .  .r.  JS(m;n  ,r)  if  the 

i  1  k  ~ 
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N  cells  are  distinguishable,  and  is  equal  to  S(m;n,r)  if  only  cells 
belonging  to  different  groups  are  distinguishable  (and  cells  in  the  same 
group  are  alike). 

It  is  easily  concluded  from  Proposition  4.1,  or  from  (*4.3)  -  (4.5); 
that  the  numbers  S(m;n,r)  are  non-negative  integers. 


5*  RATIOS  OF  MULTI  PARAMETER  STIRLING  ;jU.VLERS_OFTHE  _SJCCOKn  JOND. 

We  define,  as  ratio  of  nultiparameter  Stirling  numbers  of  the  second 

kind  with  respect  to  argument  m,  the  function 
S(mj-l;n,r) 

R2(,n;n,r)  =  H-D 

Working  as  for  Proposition  3.1fve  obtain 


Proposition  5.1.  A  recurrence  relation  for  the  ratio  B- ( r. ;  n ,  r ) ,  inde :  r . 5 •-  nt 

of  the  multiparameter  Stirling  numbers  of  the  second  kind,  is  given  by 

,  (  m  r.  I  n .  r.+l-r'n 

k  Vr.-Jly  i  3  -  - 

.  L  — J-(VV) - JT 

.  3=1  -  -  r^  i-l  J  J 


-  (5.2) 


*s. 

|T  R  (m-i,n,r) 


for  n  >  1  and  m  >  r'n  with 


R2(m,l ,r)  -  k 


(5.3) 


k  n.  k 

R?(r,n;n,r) =  (r'n+l )  n  (r.!)  1  T 


n  (r  !)  2  T  - - - 1 - -r - 

i-i  3  <1-|  r.,-1  k  n. 

‘  (r  !)  ‘‘  (mi)!  n  (r.l)  1 

0  J  i=l  1 


(5.B) 


The  special  case  k  =  1  yields 

Piojosition  5.2.  A  recurrence  relation  for  the  ratio  R^Cm.n.r), 


independent  of  the  generalized  Stirling  nuo.Lcrs  of  the  second  hind,  is 


given  by 


R^(m,D ,r)-n  = 


for  n  >  1  and  n,  >  rn  ,  with 


(rn)r 


ni+l-rn 


n 


R2(rr.-r+l-i  ,n-l  ,r) 


tn-rn 

I  1  R,(m-i,n,r) 

i=l 


{5.5) 


R2(in,l,r)  ~  1 


(5.6) 


and 


(5.7) 


R^(rn ,n ,r)  =  n(rn+l)/(r+l)  . 

Also  for  k  =  2,  r  1  we  obtain 

Proposition  5.3.  A  recurrence  relation  for  the  ratio  R^^.n),  independen 

of  the  usual  Stirling  nurbers  of  the  second  kind,  is  given  by 
r+l-n 

JT  R2(r.-i,n-l) 

R^  (rr.  ,n)-n  =  - -  (6.6) 

T~f  R„(ir,-i,r) 
i=l 

for  n  >  1  and  m  >  n  ,  with 


R2(t7.,1)  =  1 


(5.9) 


and 


R^Xn.n)  =  n (n+1 )/2  . 


(5.10) 

Proposition  5.4.  An  alterative  recurrence  relation  for  the  ratio  R2(ir.,n,r) 
is  given  by 


_  .  .  m  i-  . 

R0(ir.,n,r)-r.  - - r~ -  c 

2  *  ir.-r+l  R^r.-l  ,r.  ,r) 


£p2  (r.-l  ,r.,r)-nj  R2(jn-r,n-l  ,r) 


(5.11) 


for  n  >  1  and  m  >  rn. 

Applying  Proposition  5.U  with  r  =  1  gives 
Proposition  5.5.  An  alternative  recurrence  relation  for  the  ratio  R  (r.,n), 
is  given  by 

jR0(n-l.,n)-n  |  R0(it;-1  ,n-l ) 

(6.1?) 


|R»(n-l.,n)-n  I  R  (m-l,n-l) 
R„(rr,,n)-n  =  -- -  2 


R2(r,-l,n) 

f  or  n  >  1  and  tt.  >  n . 

The  last  relation  was  also  derived  by  Berg  (l?75)* 
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6.  MULTI  PARAMETER  C- NUMBERS , 


i  The  r.ult  iparan.eter  C-numbers,  C(rr.;n  ,s  ,r) ,  are  defined  by  their  egf 
•  s_  n  k  ~  r  s  ri"3  “I  ni 

=  L,  c(r**.r:'*!'i)  h=  rr  -j-r  l]1)  t3J  (6-^ 

-  m=r'n  2=  1  1  “  j^O  J 

where  the  s.  /  0  i  =  1 . k  are  any  real  numbers. 

Taking  k  =  1  gives  the  generalized  C-numbers,  Charalambi des  (!97^&), 
and  k  =  1 ,  =  3  defines  the  simple  C-numbers ,  Cacoullos  and  Charalambides 

(irMb),  Charalarbides  (1977). 

Tn«  following  properties  of  C(n,n,s,r)  are  easily  verified. 


a)  They  have  the  representation. 


£  ft  ft  Q,) 

where  the  summation  extends  over  all  ordered  h'- tuples  (K=n^+  •  •  •  +ny  )  of 


(6.2) 


inters  >:jJ  outi  r  Tying 


i  n  • 

k  i 


*ii  -ri  5  al»****k  and  II  Jl*  ij*  m  * 

t)  They  satisfy  the  following  recurrence  relations, 

k  i  . 


C(rtl;n,s,r)  =  (s  ’n-m)  C(m;n  ,s  ,r)  ♦  II  (ft-  J  ( s ;  )  r  C(n.-r .  h  l;n-e,  ,s  ,r) ,  (6.3) 

— 2  j  r^  i  -v  -vi  -w.  -w 

n.  (m)^  /  Si  i  j 

C(r.;n,s,r*e.)  =  (-1)3 - vft—  (  _  /  C(m-jr.;n-je  ,s,r),  (6. A) 

- 1  -i_n  \r,-  '  2  -  - 


C(r.;n,s,rte.)  =  (-1) 

~  ~  j=0 

with  initial  conditions 


C(C;C,s,r)  =  j,  C(n.;n,s,r)  =  0  when  it.  <  r'r..  They  are  obtained  from  the 

di f f er  ence-di f  f er ent i al  equation 

.  k  (s.)r.  r.-l 

(Itt)  cpn(t;s,r)  =  s‘n  (pr(t;s,r)+  J~  t  1  <pn_e  (t;s,r).  (6.6) 

-  “  ~  ""  ~  i=l  '  i  ~  -i 

The  representation  (6.2)  leads  to  the  following  int erpret at  ion  of  the 
C(n  Jn  ,s  ,r  )-rruml  ers  in  the  framework  of  coupon-collecting  problems. 

Consider  an  urn  containing  k  groups  (sets)  of  distinguishable  balls; 
th<:  ith  group  consists  of  s^.ru  balls  and  is  divided  into  equal  subgroups 
(su!  rets)  of  s^.  balls  each  bearing  the  numbers  l,...,n^;  moreover,  suppose 
that  the  halls,  of  the  k  groups  are  distinguished  by  different  colours  so 
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that  each  hall  in  the  urn  is  dist inguished  by  its  colour  and  number-  Now  it 
is  easily  seen  iron.  (0.2)  that  ^ 

Proposition  6.1.  The  number  of  ways  of  selecting  m  balls  out  of  an  urn 
k 

with  s'n  s.n.  distinguishable  balls,  divided  into  k  groups  by  colour 

i=l 

and  number  as  above  into  rm  subsets  of  size  s,.  within  the  i-th  subgroup,  so 
that  each  number  T,...,n^  of  the  i-th  subgroup  (colour)  appears  at  least 
times  is  equal  to 


Uj  .  ...  . 

- ; — - —  C(m;n,s,r). 


(0.0) 


Here  it  was  assumed  that  s^  is  a  positive  integer.  If  is  a  negative 

integer,  say,  =  -s*  ,  then 

v  X) «  ■,  n .  4  <»  ri .  .  * 

y  2  /  5  .  |  K  2.  /  “£  .  \  X  3  X../S.+X..-l\ 

nr (x*  -tt  rt  -rr  it  <-»i5(  y5  )  «■.*> 

i=l  j=l  1  ijJ  i=l  5=1  i=l  j=l  1  xij  / 

and  from.  (0.2)  it  can  be  concluded  that  the  sign  of  C(m,n,s,r)  is  the  sane 

as  (-l)n'.  Furthermore,  we  ir.ay  deduce 

Proposition  6.2.  The  number  of  ways  of  distributing  m  (m  >  r’n)  non- 
distinguishable  balls  into  s*  n  cells,  divided  into  k  groups  of  cells  with 
Siri  ce^s  t^ie  group  and  ik  subgroups  each  of  s^  cells  ir.  t).e  i-th 

group,  so  that  each  subgroup  of  the  i-th  group  contains  at  least  balls  is 
equal  to 


- - —  j  C(ir.;n,-s»  ,r)J  . 


(0.6) 


As  an  indication  of  the  applicability  of  the  nultiparameter  C-nunbers 
in  occupancy  problems,  we  refer  to  a  problem  posed  by  Scbel  et.  al .  ( ly  ) ,  \  .  ‘,2. 
Signless  multiparanreter  C-nunbers.  From,  the  basic  recurrence  relation  (6.3) 
or  from  the  last  two  propositions,  we  conclude  that 

(i)  for  s^  >  0  integer,  the  numbers  C(m-,n,s,r)  are  non-negative  integer's;  they 
are  positive  for  r’n  <  it.  <  s’n;  otherwise  sero. 

(ii)  for  s^  <  0  integer, the  numbers  C(m.;n,r,,r)  are  integers  having  the  sign 

of  (-if. 
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Thus,  as  in  the  case  cl  the  Stirling  tiuriers  of  the  first  kind,  Riordan  (tc,  >), 
the  positive  numbers  ^ 

|c(m;n,-s*,r)|  =  f-l)™  C(m;n  ,-s*  ,r)  (6.9) 

**  v  *v  <kv 

will  be  called  signless  multiparameter  C-numbers. 

It  can  be  easily  verified  that 

Proposition  6.3.  The  egf  of  the  signless  mult iparameter  C-numhers  |  C(iv,n,-s  ,r ) |  , 

s.  >0  i  =  1 . k  is  given  by 

1  k  r  -  ri”3  -s  1  ni 

wt-(z;-s,r)  -  n  r?  Si-  £  ^  (  i  )  ,j  J  \ 

n -  i=1  L  5^0  J 

Ranark.  It  should  be  observed  that  this  is  exactly  the  egf  required  for  the 

treatment  of  the  rr.vue  problem  in  the  negative  binomial  case  when  the  prob¬ 
ability  function  of  the  i-th  sample  is 

/S.4>-..-l\  X..  S.  X../-S.  \  X..  s 

P(x  =  *..)  =  —  rj - T  (  1  13  e  2=  (-D  1J  (  )e  1;,(i-e) 

g(6,r.)  \  Xij  ’  \  x^/ 

(6.11) 


g(e,r.)  =  (i-e)  1  ^  f  *a)  ej  ,  i  =  i,...,k. 

1  o  3 


(6.12) 


7.  RATIOS  OF  MULTI PARAMETER  C-NUMBERS . 

Ve  define,  as  ratio  of  mult iparameter  C-numbers  with  respect  to  argument 


m,  the  function 


C(m+1 ;n ,s ,r) 


R„  (m;n  ,s  ,r)  =  — : - r 

3  -  *V.  C(m;n,s,r) 


(7.1) 


Proposition  7.1.  A  recurrence  relation  for  the  ratio  K^(m,n,s,r),  independent 
of  the  multiparameter  C-numbers,  is  given  by 


k  (  (s . )  n .  nwl-r’n 

/—  (r'n)  /s;  1  .  ' 


W  (b) 

R3U,-,tl,s,rh-„-l_’T, - V-V'g 


R  (m-r. tl-i ,n-e . ,s,r) 

"  D  —  —  J  **  ~ 

- -  (7.2) 


IT  K,(m-i  ,n,s,r) 
i=l  3 


for  n  >  1  and  m  >  r’ri,  with 


R.  (n.,1  ,s  ,r)  =  r-n. 

v  (r’i.i+1)  k  r  nr3  k  r  n 

- -,7-  h\<<r  T)(r‘  }  ‘ 

n  irJ)  J  ^  •'  J.;j 


(7.3) 


i  - 1  ri 


(7.4) 
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Proposition  7.2.  A  recurrence  relation  for  the  ratio  R3(jr.,n  ,s  ,r) , 

independent  of  the  generalized  C-numbers  (cane  k=l),  is  given  by 

l <*>„ »  "'Vr™ 

— -  -  7T  R  (m-r+l-i.n-l.s.r) 


»s  *r)+ni-sn= 


for  n  >  1  and  n.  >  rn,  with 
R3(m,l ,s ,r)  =  s-m 

and 


jii-ra 

tt 

i=l 


TT  RgCm-i  .n.s.r) 


R  (rn,n,s,r)  =  n(rn+l)(s-r)/(r+l) 

o 

Proposistion  7.3.  A  recurrence  relation  for  the  ratio  R3(ir.,n,s) 
independent  of  the  usual  C-numbers  (case  r=l),  is  given  by 


R3(ir.,n,s)+r-sn  = 


ir,+  l-n 

TT  R,(»T.-i,n-l,s) 
i=i  3 


ir.-n 

7T  RAir.-itn,s) 
i=l  3 


for  n  >  1  and  m  >  n,  with 
R3(tt.,1  ,s)  =  s-m 


and 


R3(n,n,s)  =  (s-l)n(n+l)/2  . 


Proposition  7.4.  An  alternative  recurrence  relation  for  the  ratio 

^RgCm-l  ,n  ,s  ,r)+m-sn-lj  R3(r.-r,n-l  ,s  ,r) 


R3(m,n,s,r),  is  given  by 


m 


R3(ir,,n,s  ,r)tr.-sn  =  Ir_r+1 


R3(m-l,n,s,r) 


for  n  >  1  and  m  >  rn. 

Proposition  7.5.  An  alternative  recurrence  relation  for  the  ratio 


R3(ir.,n,s),  is  given  by 
R~(ni,n,s)-»i:.-sn  = 

w 


^R3(m-1  ,n  ,s  )tir.-sn-lj  Rg(  it,- 3  ,n-l  ,s) 


R3(r.-1  ,n,s) 


for  n  >  1  and  it.  >  n. 


(7.5) 

(7.6) 

(7.7) 

(7.8) 

(7.9) 

7.10) 

(7.11) 

.1?) 
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8.  RELATIONS  BETWEEN  THE  STIRLING  AND  C-NUMBERS. 

It  was  observed  in  Cacoullos  and  Chara lamb ides  (l gf^),  that 
ljrn  s  m  C(m,n,s)  =  S(m,n) 


(8.1) 


that  is,  the  C-numbers  can  be  approximated  by  the  Stirling  numbers  of  the 
second  bind  for  large  s,  a  fact  which  reflects  the  corresponding  well  known 
convergence  of  the  binomial  to  the  Poisson  (s  -*  »  ,  p  -*0,  ice. 6  =  p/q  ■*  0 
and  hence  sp  or  s6  converges  to  the  Poisson  parameter  X  ).  The  above 
property  extends  to  the  case  of  multiparameter  Stirling  numbers  of  the  second 


kind  and  mult (parameter  C-numbers,  namely, 

lim,  s.m  C(m;n,s,r)  =  S(m;n,r),  i  =  l,...,k. 

S  .  00  3.  «  ^  iw 


(8.2) 


This  can  easily  be  verified  by  using  the  corresponding  representation 


(4.2)  and  (6.2)  of  these  numbers  and  noting  that 


lim  s  k  (  ^  )  =  1/kJ . 


(8.3) 


A  relation  between  the  signless  multiparamieter  Stirling  numbers  of  the 
first  kind  and  the  multiparameter  C-numbers  reflects  the  limiting  relation 


between  the  negative  binomial  and  the  logarithmic  series  distributions: 


lim  s.N  |c(m;n,-s,r)|  =  |s(m',E,r)|  ,  N  =n  +...+n 

n  1  ^  ^  1  r 


(8.4) 


s  0 


1This  can  be  seen,  eg,  by  showing  that  the  egf  of  the  s.  |c(m,n,-s ,r) | - 
numbers  converges  to  the  egf  of  the  |s(m;n,r) | -numbers,  that  is, 

V.io  ~ r  JT  nTT  D1-'’  Si‘  t  ‘-n5  (  Js)  JT-hf-ioEa-t)-!:  £) 

i  s.  1=1  i  j=0  J  a=l  a  u  i=l  J 


Tor  this  note  that 


1  <-»3  c;)  <j 


(8. 6") 


s  -»  0 
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APPENDIX  B 

CHARACTERIZATIONS  OF  DISCRETE  DISTRIBUTIONS  BY  A  CONDITIONAL 
DISTRIBUTION  AND  A  REGRESSION  FUNCTION* 

T.  Cacoullos  and  H.  Papageorgiou 
Statistical  Unit,  University  of  Athens 


Key  Words  and  Phrases:  Characterizations, discrete  mixtures,  identifiability, 

regression. 


Abstract 

The  bivariate  distribution  of  (X,Y),  where  X  and  Y  are  non-negative 
integer- valued  random  variables,  is  characterized  by  the  conditional  distribution 
of  Y  given  X  and  a  consistent  regression  function  of  X  on  Y.  This  is 
achieved  when  the  conditional  distribution  is  one  of  the  distributions :a)  binomial, 
Poisson,  Pascal  or  b)  a  right  translation  of  these.  In  a)  the  conditional 
distribution  of  Y  is  an  x-fold  convolution  of  another  random  variable  independent 
of  X  so  that  Y  is  a  generalized  distribution.  A  main  feature  of  these 
characterizations  is  that  their  proof  does  not  depend  on  the  specific  form  of  the 
regression  function.  Moreover, it  is  shown  that  the  characterizations  hold  if  the 
regression  function  is  replaced  by  any  higher-order  conditional  moment.  It  is  also 
indicated  how  these  results  can  be  used  for  goodness-of-fit  purposes. 


*  A  shortened  version  of  this  paper  is  appearing  in  the  Ann.  Inst.  Statist.  Math. 


1 .  INTRODUCTION 


s 

Here  we  are  concerned  with  characterizing  the  distribution  of  non¬ 
negative  integer-valued  random  variables  (r.v.)  X  and  Y  in  terms  of 
the  conditional  distribution  of  Y  given  X  and  the  regression  function 
E(X|Y)  of  X  on  Y.  Several  papers  have  appeared  in  this  direction. 

Kbrwar  (1975)  considered  a  conditional  binomial  or  Pascal  distribution 
combined  with  linear  regression  and  characterized  the  Poisson,  binomial 
and  negative  binomial  distributions  in  the  former  case  and  the  geometric 
in  the  latter  case.  Dahiya  and  Korwar  (1977)  extended  these  characterizations 
to  bivariate  X  and  Y  under  conditional  distributions  which  are  independent 
binomials  or  Pascal  and  linear  regression.  Khatri  (1978a) ,( 1973b) .using  a 
slightly  more  general  approach  gave  similar  results  for  the  multivariate 
case.  A  case  of  non-linear  regression  was  treated  by  Xekalaki  (1980)  in 
characterizing  the  bivariate  Poisson  distribution. 

In  this  paper  more  general  and  unifying  results  are  obtained  by  by¬ 
passing  the  unnecessary  details  involved  in  obtaining  a  specific 
characterization  under  a  specific  regression  function.  This  is  achieved 
by  appealing  to  the  unicity  of  a  solution  of  a  first-order  difference 
equation.  Specifically,  it  is  shown  that  certain  conditional  distributions 
along  with  the  regression  functions  determine  uniquely  the  distributions 
of  X  and  Y,  hence  also  of  (X,Y),  This  fact  can  be  used  to  generate  a 
wide  spectrum  of  distributions  characterized  under  these  conditions.  For 
example,  given  a  conditional  binomial  distribution  (of  Y  on  X),  we  may 
choose  an  arbitrary  given  distribution  for  X,  which  in  turn  gives  a  specific 
regression  function;  thus  (see  Section  ^ ) ,  in  addition  to  the  distributions 
mentioned  earlier,  characterizations  were  obtained  for  the  logarithmic 
distribution  and  several  generalized  distributions,  e.g.  ,  Neyr.an,  Poisson 

;  binomial,  binomial  Poisson,  logarithmic  binomial  etc. 

f 

1 

\ 

I 
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Finally,  it  is  shown  that  the  distribution  of  (X^f)  can 
be  similarly  characterized  if  we  replace  E(X|Y)  by  the  kth  conditional 
moment  E[x^|V]  .  for  some  value  of  k.  It  should  also  be  added  that  the 
present  approach  goes  through  in  the  multivariate  case,  which  will  be 
treated  in  a  subsequent  paper. 


*♦ 


2.  SOME  PRELIMINARIES 

We  shall  make  use  of  the  following  easily  verified  combinatorial 
identities : 

■  *y»UK)  • 

*\yY  Vy)(y“0  ,Vl(y*k-l)t'"tA0(y)(y)  • 

*k(£)=  • 

»CT>  'y*"G7i)  <T)  • 


(2.1) 

(2.1) 

(2.2) 

(2.2) 

(2.3) 

(2.3) 


where  the  A^(y)  ,B.  (y),C..(y)  denote  polynomials  in  y  of  degree  k.  In 
m  and  r  are  integers  and  r>0. 

Moreover,  a  great  role  in  the  sequel  is  played  by 

Theorem  2.1  (Goldberg  1958,  p.61).  The  linear  difference  equation 
of  order  n 


f0(k)^n4fl(k)ykUrl+",+fn-l(k)yk+l+fn(k)>\  =  E(k> 

over  a  set  of  consecutive  integer  values  of  k  has  one,  and  only  one, 
solution  for  which  valuer,  at  n  consecutive  k-valuos  .ire  arbitrarily 


prescribed . 


3.  THE  MAIN  CHARACTERIZATI(^N  THEOREMS 
Certain  forms  of  the  conditional  distribution  p(y|x)  of  Y 
given  X=x  together  with  the  regression  function  m(y)  =  e[x|y=£]  of 
X  on  Y  determine  the  distributions  of  X  and  Y.  It  is  to  be  observed 
from  the  outset  that  in  the  main  case  considered  here,  namely,  when 
p(y|x)  is  an  x-fold  convolution  of  a  non-negative  integer-valued  r.v. 
with  probability  generating  function  (p.g.f.)  h^Cv),  it  is  sufficient 
to  determine  the  distribution  of  either  X  or  Y,  This  is  so  because  the 
p.g.f. * s  g(u)  of  X,  h(v)  of  Y  and  G(u,v)  of  (X,Y)  Cacoullos  and 
Papageorgiou  (l98lb)  satisfy  (3.1)  (see  also  (3.6a)-(3-10a) ; equivalently, 
because  of  the  identif lability  of  the  mixtures  defined  by  ( 3 .6 ) — ( 3.10) ,  with 
mixing  variable  (parameter)  x  (see  e.g.,  Teicher  1961) . 

h( v)  =  g(hQ(v) ) ,  G(u,v)  =  g(uhQ(v))  (3.1) 

Here,  we  find  first  the  probability  function  p(y)  of  Y. 

A  similar  remark  applies  to  the  case  in  which  p(y|x)  determines  a 

shift  to  the  right  by  x  of  a  non-negative  integer-valued  r.v.  with  p.g.f. 

h  (v),  so  that  denoting  by  hx(v)  the  conditional  p.g.f.  of  Y  given  X=x, 

h^(v)  =  vXh  (v),  h(v)  =  h  (v)g(v).  (3.2) 

A  characteristic  of  the  following  characterizations  is  that  their 
derivation  is  independent  of  the  actual  form  of  the  regression  function 
m(y);  this  was  confined  to  be  linear  in  the  relevant  statistical  literature. 
The  simplification  is  achieved  by  erploying  a  basic  theorem  about  the 
uniqueness  of  a  solution  of  linear  difference  equation  (Theorem  2.1).  The 
difference  equation  governing  p(y)  in  each  case  is  obtained  by  using 
Lemma  3._1_.  Let  p(y|x)  be  such  that  there  exist  p.g.f.'s  h^,  h*  independe 
of  x  and  functions  cj  of  y  (only)  so  that 

h  (v)  =  [h  (v)lX  or  h  (v)  =  vXh  (v)  for  x=0,l,?,...  (3.3) 

x  L  0  -J  y. 


xp(ylx)  =  c1(y)p(y+llx)+c0(y)p(y|x)+c_l(y)p(y-l jx) 


(3.«0 
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Then  p(y|x)  and  m(y)  characterize  the  distribution  of  Y;  hence  of 
X  and  (X,Y). 

Proof .  We  have 

m(y)  =  l  xP  [X=x  |  Y=y]  =  \  *PMj  P[X=x^] 

x=0  x=0  pvy' 

which  by  (3.4)  can  be  written  as 

tn(y)  =  [c1(y)p(y+l)+c0(y)p(y)+c_1(y)p(y-l)]/p(y)  (3.5) 

Since  this  is  a  linear  second-order  difference  equation  in  p(y)  the 
assertion  follows  by  Theorem  2.1  and  in  virtue  of  (3.1)  and  (3.2). 

Note:  In  the  applications  of  the  Lemma,  either  c^(y)  or  c  ^(y)  is 
zero  so  that  (3.5)  reduces  to  a  first-order  difference  equation.  Now 
we  state  the  main  theorems. 

Theorem  3.1.  Let  the  conditional  distribution  p(y|x)  be  one  of 
the  distributions  (3.6)-(3.10) .  Let  m(y)  =  E[x|Y=y]]  be  an  arbitrary 
function  of  y  consistent  with  p(yjx).  Then  p(y|x)  and  m(y)  together 


determine  the  distributions  of  X,Y  and  (X,Y). 

p(y|x)  =  QpV  Y  ,  y=0,...,x,  x=0,l , . . .  (q=l-p)  (3.6) 

p(y|x)  =  (x_i)pxqy_x  »  y=x,x+l,...,  x=l,2,...  (3.7) 

p(y|x)  =  (X+y_1)(|f  (1-  |)X.  y=°»1 .  X=l,...  (P=Q-1)  (3.8) 

p(y|x)  =  e  *x  ,  y=0,l,...,  x=0,l,...  (3.9) 

p(y|x)  =  QX^p^qnX  y  y=0,l, —  ,nx  x=0,l,...  (3.10) 


Proof .  By  Lemma  3.1,  it  is  enough  to  show  that  for  some  p.g.f.  h0(v) 

hx(v)  =  [vv)r 


and,  moreover,  that  p(y|x)  satifies  (3.4). 
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Indeed,  the  hQ(v)  corresponding  to  the  distributions 

\ 

(3.6)-(3.10) 

the  respective  p.g.f.'s  of  X  and  Tf  (see  (3.1))  are 

as  follows: 

hQ(v)  =  pv+q,  g(v)  =  h(^^) 

(3.6a) 

h0<v)  =  •  B<v)  =  h<p,qv) 

(3.7a) 

h0(v)  =  (Q-Pv)  ,g(v)  =  h(Qpv  * 

(3.8a) 

hQ(v)  =  e  ,  g(v)  -  h(—  +1) 

(3.9a) 

1  /n 

(3.10a) 

h0(v)  =  (pv+q)n ,  g(v)  =  h( — 

As  regards  (3.4),  by  using  (2.1)  for  (3.6)  and  (3.10),  (2.2)  for 
(3.7)  and  (2.3)  for  (3.8),  we  obtain  the  respective  difference  equations 
(recurrences)  of  the  first  order  : 


m(y)  =  %y+l)  ^ y  +y 

y=0,l,... 

(3.6b) 

m(y)  =  y-q(y-l)  Pp^~y 

y-1 ,2, . . . 

(3.Tb) 

m(y)  =  (y+Dp  P^fyy~  -  Y 

y-0 ,1 , . . . 

(3.8b) 

«<»>  *  *T  1 

y=o , i , . . . 

(3.9b) 

m(y>  =  3_  (y+1)  .  * 

y=o »!*••• 

(3.10b) 

Hence  the  proof  of  the  theorem  is  complete. 

The  explicit  solutions  of  the  above  difference  equations 

(3. 6b) -(3. 10b) 

are  respectively  : 

p(y)  =  P(0)(^yT[kTT  > 

y=i » 2, . . . 

(3.6c) 

p(y)  =  pUJq'’  TT  . 

y=2 , 3» . . . 

(3.7c) 

p(y)  =  P<°>(5;  ^TkTT  <’”<k>*k>. 

y=l,2, . . . 

(3.8c) 
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p(y)  =  p(0)xy^1J  £+!r  *  y=l,2, . . .  (3.9c) 

p(y)  =  P(°^(^  fl  (nm(k)-k),  y=l,2,...  (3.10c) 


where  p(0)  or  p(l)  are  determined  from  the  condition  £  p(y)=l. 

y 

Notice  that  (3.6)  is  a  special  case  of  (3.10)  with  n=l.  Also 
that  (3.7)  and  (3.8)  correspond  to  the  two  versions  of  a  Pascal  r.v., 
which  denotes  the  number  of  failures  in  (3.8)  and  the  number  of  trials 
up  to  the  xth  success  in  (3.7). 

Next  we  consider  the  case  described  by  (3.2).  We  prove 
Theorem  3 . 2  Let  p(y|x)  be  one  of  the  distributions  (3 . ll)-(3. 13) . 
Let  m(y)  be  consistent  with  p(y|x).  Then  p(y|x)  and  m(y)  together 
determine  the  distributions  of  X,Y  and  (X,Y). 


P(ylx>  =(v"x)pyX  9n  Y+X  xsy^x+n. 


p(y|x)  =  (N+^ _1)(|)  (1"  y>x*o,  (P=Q-i) 


y-x 


.y-x 


y£x,  x=0,l,... 
N 


(3.11) 

(3.12) 

(3.13) 


Proof .  Apply  Lemma  3.1 

h  (v) 
h  (v) 
h*(v) 


i  (3.11)-(3.13)  with  h* 

(pv+q)n 

X(v-l) 

e 

(Q-Pv)“N 


given,  respectively,  by 

(3.11a) 

(3.12a) 

(3.13a) 


Now  using  the  identity 


)=  y(  n  ^-(n-y+1)  (  n  ."'j-xf  n  . 
Vy-x/  }  \y-x/  J  \y-x-l/  W-x-1 


we  can  write  m(y)  for  (3.11)  in  the  form 


-  B8  - 


For  (3.12)  we  easily  find 

m(y)  =  y~x  2fer 

and  for  (3.13),  making  use  of  the  identity 


*n:r)= » nrv^Trhnr) , 


(3.11b) 


(3.12b) 


we  have 


m(y) 


y-(Nty-l) 


P  p(y-l) 

Q  p(y) 


+ 


p  p(y-D 
Q  p(y) 


m(y-l) 


(3.13b) 


Here  again  p(y)  is  obtained  as  the  solution  of  the  corresponding  first- 
order  difference  equation.  In  point  of  fact,  we  have 


/£>y  -tV  n-k+l+m(k-l) 


»<*>  ■  in 


p(y)  =  p(0)  jy  ,|  , 


k  =1 

y  y 


m(k-l)-(N+k-l) 


.  p(°)(f  j  FT 


where  p(o)is  also  determined  since  \  p(y)=l. 

y=0 


(3-llc) 

(3.12c) 

(3.13c) 


Finally,  we  give  another  set  of  characterizations  under  the  assumptions 
of  Theorem  3.1,  i.e.  (3.6)-(3.10) ,  replacing  m(y)  =  E[x|Y=y]  by  a 
higher  condition*1  moment  m^(y)  =  E[x^|Y=y]  for  a  k>l. 

Theorem  3.3.  For  a  k>l,  let  n^(y)  be  given  function  of  y  which 
is  consistent  with  one  of  the  conditional  distributions  (3.6)— (3.10). 

Then  p(y|x)  and  m^,(y)  together  determine  p(y)  up  to  k-1  arbitrary 
probabilities,  say,  p(l ) , . . . ,p(k-l) . 

Proof .  By  using  the  generalized  identities  ( 2. 3 ) '-( 2 . 3)  '  and  working  as 
in  Theorem  3.1,  we  obtain  a  k-th  order  difference  equation  in  p(y). 


-  By  - 


Hence  the  ass^r-tion. 

Note .  The  undetermined  p(l)  , . . . ,p(k-l)  may  be  obtained  by  using  e:g.  as 
initial  conditions  the  values  of  m(y)  for  y=l,...,k-l,  since  we  always 

CO 

have  the  condition  \  p(y)=l,  thus  determining  p(0). 

y=0 

4.  SOME  APPLICATIONS -COROLLARIES  OF  THEOREM  3.1 
It  was  already  pointed  out  that  in  characterizing  the  distribution 
of  X,Y  and  (X,Y)  in  terms  of  the  conditional  distribution  p(y|x)  of 
Y  given  X=x  and  m(y)=ii[x | Y=y] ,  the  regression  function  m(y)  must 
originate  from  some  non-negative  integer-valued  r.v.  X  and  be  consistent 
with  .p(y|x).  This  is  the  spirit  of  the  following  characterizations- 
applications  of  Theorem  3.1.  In  the  applications  of  Theorem  3.1,  it  should 
be  observed  that,  for  a  specific  r.v.  X,Y  has  the  generalized  (compound) 
distribution  of  X  by  another  r.v.  Z^  (denoted  by  XvZ^),  since  Y  has 
the  representation 


Y  =  Z1+<-.+Zx 

where  the  Z^  are  i.i.d.,with  p.g.f,  hQ(v),  independent  of  X. 

The  following  characterizations,  in  addition  to  the  usual  discrete 
distributions,  cover  some  more  involved  generalized  discrete  distributions, 
of  which  the  more  interesting  ones  are  presented  here. 

Proposition  4.1.  Suppose  (3.6)  holds.  Then  we  have: 

(a)  For  some  X>0 

m(y)  =  y+Xq  ,  y=0,l,2,..., 

iff  X  is  Poisson  (X);  then  Y  is  also  Poisson  (Xp)  and  (X,Y)  is  a 
bivariate  Poisson-Bernoulli  model,  studied  by  Lei  ter  and  Hamdan  (1973)  and 
Cacoullos  and  Papageorgiou  (1980). 

(b)  For  some  0<p'=l-q'<l  and  some  integer  n>0 


m(y)  = 


4  -r- 

+p  q 


y+n 

q  +p  q 


0^y<n . 


-  BIO  - 


iff  X  is  binomial  (n,p');  then  Y  is  also  a  binomial  (n,pp')  and 
(X,Y)  is  a  special  case  of  a  bivariate  binomial  distribution. 


(c)  For  some  N>0,  P=Q-1>0 

,  v  Q  ^  NPq  n  . 

m(y)  =  0^y  +  Q^Pq  Y'0*1’-  •  • 

iff  X  is  negative  binomial  NB(N,P);  then  Y  is  also  NB(N,Pp)  and 
(X,Y)  is  a  bivariate  negative  binomial-Bernoulli  model,  studied  by 
Cacoullos  and  Papageorgiou  (198la). 

(d)  For  some  O<0<1 

=  y  y=1»2**** 

-  1  0q 

m(0)  “  logd-eqV  l-6q 

iff  X  is  logarithmic  (6);  then  Y  follows  a  modified  logarithmic  (6,6"), 
with  6  the  probability  of  Y=0  and  6'  its  ordinary  parameter,  where 


•■pWjS&SJ1  » 


(X,Y)  is  a  bivari-te  logarithmic-Bernoulli  model  (cf.  Cacoullos  and 
Papageorgiou ,  ( 198lb ) . 

Remark .  Cases  (a),  (b)  and  (c),  m(y)  =  ay+b,  y=0,l,...  ,  as 

considered  by  Korwar  (1975).  Case  (d),  exhibiting  linearity  only  for 
y>0,  was  missed  by  Korwar  in  view  of  the  limitations  of  treating  the 
problem  in  terms  of  a  specified  regression  function.  The  present  more 
general  approach  allows  a  variety  of  m(y),  including  those  which 
correspond  to  truncated  versions  of  X.  For  example,  if  X  is  truncated 
on  the  left  at  r,  m(y)  is  linear  only  for  y>r  under  (3.6)  and  the 
corresponding  characters  nations  hold  with  appropriate  modifications  for 
Y  and  (X,Y). 


Bll  - 


(e)  Foj?  some  positive  X  and  6 
1»<y>  .  y 

Sy(cl> 


y=o,i , — 


where  S^Ct)  is  a  Stirling  polynomial  defined  by  Charalambides  (1977), 


S  (t)  =  £  S(n,k)t  ,  c-  =  Xe  , 

k=0 


and  S(n,k)  denotes  a  Stirling  number  of  the  second  kind, 

iff  X  has  a  Neyman  distribution,  i.e.,  Poisson  (X)  vPoisson  (9); 

then  Y  has  also  a  Neyman:  Poisson  (X)  vPoisson  (0p)  and  (X,Y)_, 


with  pgf 


G(u,v)  =  exp[X{e6  Lu^+Pv^  3  _j}J  f 


is  a  special  case  of  a  bivariate  Neyman  Type  1  distribution,  Holgate  (1966) 
(f)  For  some  0<p'=l-q'<l,  X>0  and  an  integer  n>0 


'q  Cy+l,n(c?)  .  , 

C  (cj  +  y  *  y-0’1*--- 


m(y).  =  -£-4-  - 

q  *P  q  >yfn^ 


where  C  (t)  is  the  polynomial,  Charalambides  (1977), 

y 

Cv  n(t)  =  ^  C(y,k,n)tk  ,  c  =  X(q'  +  p'q)n  , 

y’  k  =  0 

and  C(y,n,k)  are  the  C- numbers  ,iff  X  is  Poisson  (X)  v  binomial  (n,p'); 
then  Y  is  also  a  Poisson  (X)  v  binomial  (n,p'p). 


It  is  worth  noting  that  (X,Y)  with  p.g.f. 


G(u,v)  =  exp[x{  fq  '+p  'u(q+pv)]n  -l)J 


is  a  special  case  of  the  bivariate  Poisson  binomial  type  I  distribution 
studied  by  Charalabmider.  and  Papageorgiou  (1981). 


-  B12 


(g)  For  some  P=Q-1>0,  A>0  and  an  N>0 


where 


=  r  <m  c^y-  y=0,1,J" 


Cy_.N(t)  =  l  C<y,k,->»tk  ,  '3  =  >(Q-Pq)',J  , 


iff  X  is  Poisson  (X)  v  negative  binomial  (N,P);  then,  Y  is  also 
Poisson  -(A)  v  negative  binomial  (N,Pp). 

(h)  For  some  0<p'=l-q'<l;  A>0  and  an  integer  n>0 

XqS  .  (ca) 

m(y)  =  s  (c~S -  +  y  y=0,l,2,... 

y,n  4 


where 


S  (t) 


=  \  (n).  S(y,k)tk,  c  =  p  *e  *P(q'+p'e  Ap)  1 


iff  X  is  binomial  (n,p*)  v  Poisson  (X);  then  Y  is  also  binomial 
(n,p*)  v  Poisson  (Xp). 

(i)  For  some  0<p '=l-q '<1 ,0<p*=l-q*<l  and  positive  integers  n,n* 


where 


m(y)  = 


c  *<C-) 

: -  +  y  ,  0<y$nn*  , 


q*+P"q  cy,n,n*(c5) 


C  ...(t)  =  l  (n).  CCy.k.n*)^ 
y,n,n-  v_n  k 


c5  =  p'(q*+p*q)n  {q '+p '(q^'+p^q)11  )  > 


iff  X  is  a  binomial  (n,p')  v  binomial  (n*,p*)»  then  Y  is  also  a  binomial 
(n,p')  v  binomial  (n*,p*p) 

(j)  For  some  P=Q-1>0,  0<p'=l-q'<l  an  integer  n>0  and  an  N>0 


C  ^  M(c,) 

,(y)  -  y-  PQ - ^-44- 

Hy'  y  T?=pq  c  ^  MfHTT 

y,n,-N  6 


y=0,l,?,. 


-  B13  - 


where 


and 


c  M(t)  =  i  (n).  Cty.k.-lOt* 
y,n,-N  k^0  k 


c6  =  p'(Q-Pq)"N{q'+p'(Q-Pq)'N)  \ 


iff  X  is  binomial  (n,p')  v  negative  binomial  (N,P);  then  Y  is 
also  a  binomial  (n,p')  v  negative  binomial  (N,Pp). 

(k)  For  some  P=Q-1>0",  A>0  and  N>0 

Aq  S  .  M(c7) 

m(y)  =  ,  y  /fv-  +  y  y=0,l,2,... 
y,-Nv 


where 


and 


S  «c.-t)  =  \  (-l)k(-N)k  S(y ,k)t* 

y’  k=0 


c7  =  Pe  *P(Q-Pe  *P)  1  , 


iff  X  is  negative  binomial  (N,P)  v  Poisson  (A);  then  Y  is  also  a 
negative  binomial  (N,P)  v  Poisson  (Ap). 

U)  For  some  0<p  '=l-q  '<1  ,P=Q-1>0 ,  N>0  and  an  integer  n>0 


m 


,  C  ..  ..  C-c_) 

M  K  y,-N,n  8 


where 


and 


C  „  (-t)  =  f  (-l)k(-N)k  C(y,k,n)t* 

y*'N»n  k=o 


cB  =  P(q '+p 'q)n{Q-P(q '+p 'q)n) 

iff  X  is  negative  binomial  (N,P)  v  binomial  (n,p');  then  Y  is 
also  a  negative  binomial  (N,P)  v  binomial  (n,p'p). 

(m)  For  some  PrQ-l>0,  P*=Q*-1>0  and  N>0,  N*>0 

„(y,  s  y  .  r-o .1  . 


-  Biir- 


-  B15  - 


iff  X  is  logarithmi^  (0)  v  binomial  (n,p');  then  Y  is  also  logarithmic 
(0)  v  binomial  (n,p'p). 

(p)  For  some  O<0<1,  P=Q-1>0  and  n>0 


C*  (c  ) 

,  S  pq  V+1,-Nvc12; 

=  y-  QT-fe  c*~  vg-  \  - 

W  q  y,-Nvc12' 


y=i»2,. 


m(0)  = 


_ E3_  c*  (c  ) 

log(l-e(Q-Pq)-N)  Q-Pq  1*-»  12  ' 


where 


■  ! 

*  k=l 


(k-1)!  CCy ,k  ,-N)t 


=  6(Q-Pq)'N{l-0(Q-Pq>  N)  , 


iff  X  is  logarithmic  (0)  v  negative  binomial  (N,P);  then  Y  is  also  a 
logarithmic  (0)  v  negative  binomial  (N,Pp)  , 

Proposition  4.2.  Suppose  (3.7)  holds.  Then  we  have: 

(a)  For  some  0<p'=l-q'<l 


m(y)  =  ?Wy+  ?Vq  y=1’2’—' 

iff  X  is  geometric  (p*);  then  Y  is  also  geometric  (p'p). 


(b)  For  some  O<0<1 


(q*6p)y  1-qy  1 
(q+6p)y-qy 


m(y)  =  y-q  -q*6p  -  -~q -  »  y=l,2,  •  •  • 


iff  X  is  logarithmic  (0);  then  Y  is  a  logarithmic  (0)  v  geometric  (p). 


Proposition  4.3.  Suppose  (3.8)  holds.  Then  we  have  : 
i  (a)  Tor  some  P=Q-1>0  and  a  X>0 

C  .  i (c. ,) 

m(y)  =  -  -l^Tc-V  ‘  y  y=0,l,2,... 
y,-l  13 

where  C  ,(t)  irr  defined  under  Proposition  ** .  1  (g)  and 

y.-i 


-  Bl6  - 


c13  =  ‘O'1  . 


iff  X  is  a  Poisson  (A);  then  Y  is  a  Poisson  (A)  v  geometric  (1/ Q)< 
(b)  For  some  0<p=l-q<l  and  an  integer  n>0 

m(y)  =  -  -y-~ — 4  -  y  ,  y=0,l,2,..., 

Cy,n,-l(c14) 

where  C  . (t)  is  defined  under  Proposition  4.1  (j)  and 

C14  =  pQ~1(q+pQ“1)“;1  , 

iff  X  is  a  binomial  (n,p);  then  Y  is  a  binomial  (n,p)  v  geometric 
(1/Q). 

Cc)  For  some  P  *=Q  '-1>0 

C 


_  'y+l,-N,-l(  c15) 

m(y)  z  -  — Q  '  ' /  -c— )  -  -  y  > 

y»-N,-iv  cis' 


y= 0,1*2,... , 


where  _N  ^(-t)  is  defined  under  Proposition  4.1  (m)  and 


ci5  =  '•-Q'1tQ--p-Q':1r1. 


iff  X  is  a  negative  binomial  (N,P*)j  then  Y  is  negative  binomial 
(N  ,P ')  v  geometric  (1/Q). 

(d)  For  some  O<0<1 

C*  ,(C-r.) 

°,<y)  -  -  -c%— ’i*  --  -  y  .  y=i. 

y,-l'  16’ 


m 


(0)  = 


-c 


16 


logd-GQ-1) 


where  C*  ^(t)  is  defined  under  Proposition  4.3  (p)  and 


ci6  = 


k 


iff  X  is  a  logarithmic  (f);  then  Y  is  a  logarithmic  (0)  v  geometric 

(1/Q). 


-  B17  - 


Proposition  4.4.  Suppose  (3.9)  holds.  Then  we  have: 
\j  (a)  For  some  X>0  and  some  0>O 


S  . (c« n  ) 
,  .  y+1  17 

m(y)  S  (c'T 

V  17 


c1?  =  6e  ,  y=0 ,1 ,2, . . . 


where  S^(t)  is  defined  under  Proposition  4.1  (e) 


iff  X  is  a  Poisson  (6);  then  Y  is  a  Poisson  (e)v  Poisson  (X)  and 
(X,Y)  is  a  bivariate  Poisson-Poisson  model  studied  by  Leiter  and 
Hamdan  (1973)  and  Cacoullos  and  Papagoorgiou  (1980) . 

\J  (b)  For  some  0<p=l-q<l  and  an  integer  n>0 

Sy+l,n(c18)  „ 

m(y)  -  /  \  ’  y-0,1,2,. ... 

Sy,n(c18) 

where  S  (t)  is  defined  under  Proposition  4.1  (h)  and 
y,n 

-X,  A  -X.-l 
clg  =  pe  (q+pe  )  , 

iff  X  is  a  binomial  (n,p);  then  Y  is  a  binomial  (n,p)  v  Poisson 
(X). 

(c)  For  some  P=Q-1>0  and  an  N>0 

S 


m 


,  \  _  "y+l,-N(  ci9  ,  _ 

(y)  -  -  )  5  y-0,1,2,..., 

y,-NV  c19 ' 


where  S  „(-t)  is  defined  under  Proposition  4.1  (k)  and 

y  ,-n 

cig  =  Pe~X(Q-pe”X)_1 

iff  X  is  a  negative  binomial  (N,p)  then  Y  is  a  negative  binomial 
(N,P)  v  Poisson  (X) . 

(d)  For  some  0<B<1  . 

.(V)  =  -’^0-  . 


m(0)  = 


'20 


loC(l-0e  X) 


-  Bl6  - 


where  S*y(t)  is  defined  under  Proposition  4.1  (n)  and 

<=20  ' 

iff  X  is  a  logarithmic  (6);  then  Y  is  a  logarithmic  (6)  v  Poisson  (X) 

Proposition  4.5  Suppose  (3.10)  holds.  Then  we  have: 

(a)  For  some  X>0 

C  (c  ) 

.  .  y+l  ,n  C21 '  v  .  n  __  ,  . 

m(y)  "  nC  (c_.)  n  *  C?1  '  Xq  ’  y-0*1’2 . 

y,n  21 

where  C  (t)  is  defined  under  Proposition  4.1  (f) 

y.n 

iff  X  is  a  univariate  Poisson  (X);  then  Y  is  a  Poisson  (X)  v  binomial 
(n,p)  and  (X,Y)  is  a  bivariate  Poisson-binomial  model  studied  by 
Cacoullos  and  Papageorgiou  (1980). 

(b)  For  some  0<p'=l-q'<l  and  an  integer  n '>0 

C  »  (c 

m(y)  =  — ^ — *-~*y — — y~  +  ^  >  0<y^<n  'n  , 

where  C  .  (t)  is  defined  under  Proposition  4.1  (i)  and 
y,r.  ,n 

c22  =  p  q  (q  +p  q  >  , 

iff  X  is  a  binomial  (n',p');  then  Y  is  a  binomial  (n^p')  v  binomial 
(n,p) 

(c)  For  some  P=Q-1>0  and  an  N>0 


C  (-c0>5) 

,  .  y+l ,-N ,n  23  y 

l(y)  '  nC  ..  +  n  ’  ' -0,1,2,... 

y,-N,n  23 


where  C  (-t)  is  defined  under  Proposition  4.1  (£)  and 

y,-N,n 


e?3  =  Pqn(Q-l'qn)'1  , 


iff  X  is  a  negative  binomial  (N,P);  then  Y  is  a  negative  binomial 
(N,P)  v  binomial  (n,p). 


(d)  For  some  O<0<1 


C*  (c  ) 

.(y):-^  24  4, 

nC'WC24>  " 


y-1 ,2,. . . , 


-C*  (c  ) 

1  ,n  c24' 

m(0)  =  - -  » 

n  log(l-6qn) 

wh_re  C*  (t)  is  defined  under  Proposition  4.1  (0)  and 

y,n 


c24  =  6qn(l-eqn)‘1> 


iff  X  is  a  univariate  logarithmic  (6);  then  Y  is  a  logarithmic 
(0)  v  binomial  (n,p). 


5.  SOME  APPLICATIONS  -  COROLLARIES  OF  THEOREM  3.2 
Proposition  5.1.  Suppose  (3.12)  holds.  Then  we  have: 

(a)  For  some  0>O 


m(y)  =  y~y  x?e  ’  y=o,i,2,...  , 

iff  X  is  a  univariate  Poisson  (0);  then  Y  is  also  a  Poisson  (0+X). 

(b)  For  some  y>0 


m(y)  =  y-y 


H*  .(6) 

y-1 


(2y)1/2  H*y(B) 


y=0,l,2,.. 


where  H*y(t)  are  the  modified  Hermite  polynomials,  Kemp  and  Kemp  ( 1  9 6 5 X 
defined  as 


[l/2y]  ylty  2d 

H*  (t)  =  l  - r 

y  j=0  (y-2j ) !j !2^ 

-1  /o 

and  B=X(2y)  . 


iff  X  is  a  doublet  Poisson  (i.e.  X=?Z  where  Z  is  a  Poisson  ( u ) ) ; 

then  Y  is  a  Hermite  distri but i on , i . e . ,  a  convolution  of  an  ordinary 

Poisson  with  a  doublet  Poisson. 

(c)  For  some  0>0  and  some  p>0 

A  ,  _jU,e,S.(iio‘0)) 

m(y)  -  y-  ■ N  -  J  _e - •  y-o,i,2,.. . 

Ay{X,0,S.(ye  )) 


-  b2  0 


where 


Ayn,6,s.(t))  =  s3,^ 


and  Sj(t)  was  defined  under  Proposition  4.1  (e), 

iff  X  is  a  Neyman  (u,0);  then  Y  is  a  "Short"  distribution  (X,vi,0), 

Kemp  (l967),i-e>»  a  convolution  of  a  Poisson  (A)  with  a  Neyman  (p,0). 

6.  SOME  STATISTICAL  APPLICATIONS  OF  THE  CHARACTERIZATIONS 
The  preceding  results,  in  addition  to  their  probabilistic  interest, 
can  be  used  in  goodness-of-fit  tests  in  a  variety  of  situations. 

For  illustration,  consider  the  case  in  which  records  (X)  of  accidents 
and  corresponding  fatal  accidents  (Y)  are  available  for  a  series  of 
periods.  Then  we  may  be  faced  with  identifying  the  distribution  of  X 
and  Y  under  the  natural  assumption  that  Y  given  X  is  hinomial. 

This  is  the  situation  described  by  (3.6).  A  possible  test,  within  the 
framework  of  these  characterizations,  is  to  look  at  the  regression 
function  m(y)  of  X  on  Y.  Thus,  if  m(y)  is  linear 


m(y)  =  a+by,  y=0,l,..., 

then  (Proposition  4.1  (a) ,(b) ,(c))  ,  a  regression  line  with  slope  b=l, 
shows  that  X  (hence  also  Y)  is  a  Poisson,  a  b<l  indicates  that  X 
(hence  also  Y)  is  a  binomial  and  a  b>l  suggests  a  negative  binomial 
for  X  and  Y.  On  the  other  hand,  a  line  m(y)  =  by  with  b>l  for 
y>l  and  an  isolated  point  at  y=0  (Proposition  4.3  (d))  indicates  a 
logarithmic  X. 

Similar  remarks  can  be  made  concerning  the  cases  of  more  complicated 
regression  functions,  which  ar  a  rule,  take  us  away  fro:,  the  simple 
classical  discrete  distributions.  This,  however,  is  beyond  the  scope 


of  the  present  investigation  and  we  shall  not  pursue  it  here  any  further. 


i  n  ► 
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